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Abstract. We present the results of a systematic, unbiased search for subpulse modulation of 187 pulsars performed with 
the Westerbork Synthesis Radio Telescope (WSRT) in the Netherlands at an observing wavelength of 21 cm. Using new ob- 
servations and archival WSRT data we have expanded the list of pulsars which show the drifting subpulse phenomenon by 
42, indicating that at least one in three pulsars exhibits this phenomenon. The real fraction of pulsars which show the drifting 
phenomenon is likely to be larger than some 55%. The majority of the analysed pulsars show subpulse modulation (170), of 
which the majority were not previously known to show subpulse modulation and 30 show clear systematic drifting. The large 
number of new drifters we have found allows us, for the first time, to do meaningful statistics on the drifting phenomenon. We 
find that the drifting phenomenon is correlated with the pulsar age such that drifting is more likely to occur in older pulsars. 
Pulsars which drift more coherently seem to be older and have a lower modulation index. There is no significant correlation 
found between P3 and other pulsar parameters (such as the pulsar age), as has been reported in the past. There is no signifi- 
cant preference of drift direction and the drift direction is not fou nd to be co rrelated with pulsar parameters. None of the four 
complexity parameters predicted by different emission models Ijenet & Gil2 0031 are shown to be inconsistent with the set of 
modulation indices of our sample of pulsars. Therefore none of the models can be ruled out based on our observations. We also 
present results on some interesting new individual sources like a pulsar which shows similar subpulse modulation in both the 
main- and interpulse and six pulsars with opposite drift senses in different components. 
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1. Introduction 

Despite the fact that explaining the emission mechanism of ra- 
dio pulsars has proved very difficult, this field has the advan- 
tage that we have very detailed knowledge about the emission 
mechanism from observations. We know from the very high 
observed brightness temperatures that the radio emission must 
be coherent, we know what kind of magnetic field strengths 
are involved and even the orientation of the magnetic axis, ro- 
tation axis and the line of sight can be derived from observa- 
tions. Furthermore if one can detect single pulses one can see 
that the pulses of some pulsars consist of subpulses and for 
some pulsars these subpulses drift in successiv e pulses in an 
organized fashion through the pulse window ( Drake & Craft 
1968t lSutton et alJl970l) . If one plots a so-called "pulse-stack", 
a plot in which successive pulses are displayed on top of one 
another, the drifting phenomenon causes the subpulses to form 
"drift bands" (an example is shown in the left panel of Fig. 0. 
This complex, but highly regular intensity modulation in time 
is known in great detail for only a small number of well studied 
pulsars. 
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Because the properties of the subpulses are most 
likely determined by the emission mechanism, we learn 
about the physics of the emission mechanism by study- 
ing them. That drifting is linked with the emission mech- 
anism is sugg ested by the fact that drifting is affe cted 
by "nulls" (e.g. iTavlor & Hu guejriinl [r97 it Ivan Leeuwen et all 
20021 Ijanssen & van Lee uwen 2004), where nulling is the phe- 
nomenon whereby the emission mechanism is switched off for 
a number of successive pulses. Another complex phenomenon 
is drift mode changes where the drift rate switches between a 
number of discrete values. For some pulsars there are observa- 
tionally determined rules describing which drift mo de changes 
are a l lowed from which drift mode (e.g. IWright & Fowled 
1981; Redman et aD l2005b . It has been found that the nulls 
of PS R B230 3+30 are confined to a particular drift mode 
jRedman et aljEoO^) . which further strengthens the link be- 
tween drifting and the emission mechanism. 

Another characteristic feature of the emission mechanism 
is that when one averages the individual pulses, the result- 
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ing pulse profile is remarkably stable over time (Helfand et al. 
Il975i) . Explaining the various shapes of the pulse profiles of 
different pulsars and their dependence on observing frequency 
has proven to be very complicated, so not surprisingly an ex- 
planation that is fully consistent with the overwhelmingly de- 
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tailed complex behavior of individual (sub)pulses, the nulling 
phe nomenon and the polarization of individual pulses (e.g. 
Edwards 2004) seems to be far away. In this paper we describe 
trends of the subpulse modulation we find for a large sample 
of pulsars. By doing this we determine observationally what 
the important physical parameters are for subpulse modulation, 
which could help formulating an emission model which is fully 
consistent with the observations. 

There are a few types of models that attempt to ex- 
plain the drifting phenomenon. The most well known 
mode l is the sparking gap model dRuderman & Suther land 
19751). whi ch has been extended by many authors (e. 



Chens & Ruderman 1980; Filippenko & Radhakr shnan 198 
Gil & Sendvkl l2000t Idil et al.1 l2003t lOiao et al.1 12004 mak- 



ing it the most developed model for explaining the drift- 
ing phenomenon. These models explain the drifting phe- 
nomenon by the generation of the radio emission via a rotat- 
ing "carousel" of discharges which circulate around the mag- 
netic axis due to an E x B drift. In the carousel model it 
is expected that all pulsars shou ld have some sort of circu- 
lation time. For PSR B0943+10 lEeshpande & Rankinll 19991 
iDeshpande & Rankinll200ll lAsge kar & Deshpande 2001) and 
possibly PSR B0834+06 ( lAsgekar & Deshpandell2005t) a ter- 
tiary subpulse modulation feature has been detected from the 
fluctuation properties and viewing geometry. This periodicity 
has been interpreted as related to the carousel modulation pe- 
riod (i.e. the circulation time P3), supporting the interpretation 
of the drifting subpulses being caused by a rotating carousel 
of sub-beams. The circulation times of these pulsars, as well as 
the more indirectly d erived circulation times of PSR B080 9+74 



(Ivan Leeuwen et alJ 120031) and PSR B0826-34 (Em 




are consistent with the sparking gap model (iGil et al 



A different geometry of the polar cap of PSR B0826-34 
is proposed by Esamdi n et alJl2005t In their interpretation the 
carousel changes drift direction, something what would be in- 
consistent with the sparking gap model. 

These models still have problems, like explaining the sub- 
pulse phase steps which are observed for some pulsars. Two 
clear examples of pulsars that show subpulse phase steps are 
PSR B 032 0+39 and PSR B0809+74 as fo und bv lEdwards etall 
J2003I) and lEdwards & Stapperd J2003bh . We find that the new 
drifter PSR B2255+58 also shows a phase step. 

Non-radial pulsations of neutron stars were orig inally pro- 
posed as the origin of the radio pulses of pulsars (Ruderman 



1968]) and later as a possible origin of the drifting sub- 
pulses (Drake~& Craftl Il 968t). R ecently this idea was re- 
vised bv Tciemens&Rosenl 1120041) . This model gives a natu- 
ral explanation for observed subpulse phase steps, nulls and 
mode changes. This model can be tested, although there are 
many complications, by exploring average beam geometries. 
Although this model can explain phase steps, it cannot explain 
the curvature of the drift bands of many pulsars (see Sect. 14.51 
for details). In this model it is also difficult to explain pul- 
sars with opposite drift senses in different components, because 
drifting is a simply a beat between the pulse period and the pul- 
sation time. Bi-drif ting is recently observed for PSR J0815+09 
(McLaugh hn etalJ l2004). In this paper we show a number of 
other pulsars with opposite drift senses in different compo- 



nents 1 . For PSR B1839-04 we observe that the two compo- 
nents have mirrored drift bands (i.e. the components drift in 
phase) like PSR J08 15+09, something we do not know for the 
other pulsars. In the sparking gap model bi-drifting can be ex- 
plained if thes e pulsars h ave both an inner annular gap and an 
inner core gap dOiao et all2004l) . 

A feedback model is proposed by Wright (2003) as a natu- 
ral mechanism for both the sometimes regular and sometimes 
chaotic appearance of subpulse patterns. In this model the outer 
magnetosphere interacts with the polar ca p and the obse rved 
dependency of conal type on pulse period (Rankin 1993a) and 
angle between the rotation and magnetic axis (Rankin 1990) 
follows naturally. 

Up to now most observational literature on the drifting phe- 
nomenon has been focused on describing individual very inter- 
esting drifting subpulse pulsars. The focus of this paper will 
not only be the individual systems, but also the properties com- 
mon to t he pu ls ars that show drif ting, an approach s tarted by 
iBackusI Jl98lh.lAshworthl i 1982b and iRankinl dl986h . In the 
work of Backus (1981) 20 pulsars were studied for their sub- 
pulse behavior at 430 MH z and 9 were observed to be drifting. 
In the work of lAshworthl Jl98 2) the single pulse properties of 
nine new drifters are described and the properties common to 
28 drifters in a sample of 52 pulsars are analysed. This sample 
consists of both their own results and a few previously pub- 
lished results. Most observations were obtained at or near 400 
MHz, but some at hi gher frequenc ies. 

In the work of Rankin ( 1986) all the, then published, sin- 
gle pulse properties are combined and described in the light of 
her empirical theory. Because understanding the drifting phe- 
nomenon is considered important for unraveling the mysteries 
of the emission mechanism of radio pulsars, we decided that it 
was time to start this more general and extensive observational 
program on the drifting phenomenon. 

The main goals of this unbiased search for pulsar subpulse 
modulation is to determine what percentage of the pulsars show 
the drifting phenomenon and to find out if these drifters share 
some physical properties. As a bonus of this observational pro- 
gram new individually interesting drifting subpulse systems are 
found. In this paper we focus on the 21 cm observations and in 
a subsequent paper we will focus on lower frequency observa- 
tions and the frequency dependence of the subpulse modulation 
properties of radio pulsars. 

The list of pulsars which show the drifting phenomenon 
is slowly expanding in time a s more sufficiently b right pul- 
sars are found by surveys ("e.g. iLewandowski et alJEo04l) . but 
we have successfully chosen a different approach to expand 
this list much more rapidly. The reason that we have found so 
many new drifting subpulse systems is twofold: we have ana- 
lyzed a large sample of pulsars of which many were not known 
to show this phenomenon, and we used a sensitive detection 
method. Previous studies of drifting subpulses often used track- 
ing of individual subpulses through time, an analysis method 
that requires a high signal-to-noise (S/N) ratio because it re- 



1 PSRs B0450+55, B1540-06, B0525+21, B1839-04, B2020+28, 
the outer components of B0329+54 and possibly PSR B0052+51. 
Also PSR B 1237+25 is a known example. 
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quires the detection of single pulses. This automatically implies 
that this kind of analysis can only be carried out on a limited 
number of pulsars. Analyzing the integrate d Two-Dimens ional 
Fluctuation Spectrum (2DFS; Edwards & Stappers 2002 ) and 
the Longitude-Resolved Fluctuation Spectrum (LRFS; Backer 
1970) allows us to detect drifting subpulses even when the 
S/N is too low to see single pulses. This method was already 
successfully used with ar chival W esterbork Synthesis Radio 
Telescope (WSRT) data bv lEdwards & StappersN2003al> to find 
drifting subpulses in millisecond pulsars. 

By using the technique described above combined with the 
high sensitivity of the WSRT we have analyzed a large sam- 
ple of 187 pulsars. An important aspect when calculating the 
statistics of drifting is that one has to be as unbiased as possi- 
ble, so we have selected our sample of pulsars based only on 
the predicted S/N in a reasonable observing time. While this 
sample is obviously still luminosity biased, it is not biased to- 
wards well-studied pulsars, pulse profile morphology or any 
particular pulsar charac teristics a s were previou s studies (e.g. 
Ashworth 1982, Backus 1981 and Rankin 1986 and references 
therein). Moreover, all the conclusions in this paper are based 
on observations at a single frequency. 

The paper is organized such that we start by explain- 
ing the technical details of the observations and data analy- 
sis. After that the details of the individual detections are de- 
scribed and in table |2] all the details of our measurements 
can be found. After the individual detections the statistics of 
the drifting phenomenon are discussed followed by the sum- 
mary and conclusions. In appendix |A] are the plots for all 
the pulsars in our source list. They can also be found in ap- 
pendix |B] but there they are ordered by appearance in the 
text. Note that the astro-ph version is missing the appendices 
due to file size restrictions. Please download appendices from 
http : //www . sci ence . uva . nl/~wltvr ede/2 1cm . pdf 

2. Observations and data analysis 

2.1. Source list 

All the analyzed observations were collected with the WSRT 
in the Netherlands. The telescope is located at a latitude of 
52°. 9 in the north, meaning that not all pulsars are visible for 
the WSRT. Only catalogued 2 pulsars with a declination (J2000) 
above -30° were included in our source list. 

This list of pulsars that are visible to the WSRT was sorted 
on the observation duration required to achieve a signal-to- 
noise (S/N) ratio of 130. Of this list we selected the first 191 
pulsars, which required observations less then half an hour in 
duration. The S/N ratio of a pulsar observation c an be predicted 
with the following equation jDewev etalll985l) 



S/N 



-^sys + ^~sky 



tjqS G / Ay fobs^-Po - w) 



(1) 



where t\q is the digitization efficiency factor, S the mean flux 
density of the pulsar, G the gain of the telescope, r sys the sys- 
tem temperature, r s k y the sky temperature, Av the bandwidth of 



the pulsar backend, f b s the observation duration, n p the num- 
ber of polarizations that are recorded, Pq the barycentric pulse 
period of the pulsar and w the width of the pulse profile. 

All observations were conducted with the 21 cm backend 
at WSRT, which has the following receiver system parameters: 
tjq = 1, G = 1.2 K/Jy, r sys = 27 K, T s ^ y = 6 K (which is 
the average of the entire sky), Av = 80 MHz and n p = 2. It is 
required that the pulsars have an integrated pulse profile with a 
predicted S/N ratio of 130, so the required observation duration 
in seconds is 



fobs ^ 



283 mJy \ 2 



w 



1400 



o - w 



(2) 
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where S 1400 is the flux of the pulsar at our observation fre- 
quency of 1400 MHz and w the FWHM of the pulse profile. 
Those pulsars lacking the necessary parameters (S 1400 and w) 
in the catalog were excluded from the sample, because in such 
cases it was not possible to evaluate f t, s . 

The sensitivity to detect drifting subpulses does not only 
depend on the S/N ratio of the observation, but also on obtain- 
ing a large number of pulses. This is because the observation 
should contain enough drift bands to be able to identify the 
drifting phenomenon. Our second requirement on the minimum 
observation length was therefore that the observations should 
contain at least one thousand pulses, so some pulsars had to be 
observed for longer than was required to get a S/N ratio of 130. 
To make sure that the statistics on the drifting phenomenon is 
not biased on pulse period, it is important to include these long 
period pulsars in the source list. 

Archival data was used if available and the sample of pul- 
sars was completed with new observations. The best WSRT 
data available was chosen, so for a number of pulsars the data 
greatly exceed the minimum S/N and the number of pulses re- 
quirement. This does not bias our sample of observations to- 
ward well-studied pulsars, because all the observations are long 
enough to provide a good chance to detect the drifting phe- 
nomenon. We have observations of all the sources except the 
millisecond pulsar B 182 1-24, because of the high time reso- 
lution required and the associated data storage problems. The 
observations of PSR B1823-13, B1834-06 and J1835-1020 
failed, and therefore are not included in this paper. 

2.2. Calculation of the pulse-stacks 

All the observations presented in this paper were made at an 
observation wavelength of 21 cm spread out over the last five 
years. The signals of all fourteen 25-meter dishes of the WSRT 
were added together by taking into account the relative time 
dela ys between them a nd processed by the PuMa pulsar back- 
end IVoute et aljEool . In order to reduce the effects of inter- 
ference, badly affected frequency channels were excluded. The 
frequency channels were then added together in an offline pro- 
cedure after dedispersing them by using previously published 
dispersion measures. 

To study the single pulse behavior of pulsars one usually 
converts the one-dimensional de-dispersed time series into a 
two-dimensional pulse longitude versus pulse number array 
(pulse-stack). An example is shown in the left panel of Fig. 
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Fig. 1. The left panel shows a pulse-stack of one hundred successive pulses of PSR B1819-22. Two successive drift bands are 
vertically separated by P3 and horizontally by Pi. The products of our analysis are shown for three pulsars. The top panel shows 
the integrated pulse profile (solid line), the longitude-resolved modulation index (solid line with error bars) and the longitude- 
resolved standard deviation (open circles). Below this panel the LRFS is shown with on its horizontal axis the pulse longitude 
in degrees, which is also the scale for the abscissa of the plot above. Below the LRFS the 2DFS is plotted and the power in the 
2DFS is vertically integrated between the dashed lines, producing the bottom plots. Both the LRFS and 2DFS are horizontally 
integrated, producing the side-panels of the spectra. See the main text for further details about the plots. 



H where one hundred successive pulses are plotted on top of 
one other. The pulse number is plotted vertically and the time 
within the pulses (i.e. the pulse longitude) horizontally. The 
off -pulse region is used to remove the baseline from the pulsar 
signal, making the average noise level zero. 

To correct for the pulse longitude shift of successive pulses 
the TEMPO software package 3 was used. Because the pulse 
period (P ) of the pulsar is not exactly equal to an integer num- 
ber of time sample intervals, each pulse (as it appears in the 
binned sequence) is effectively shifted by a consta nt amount 
modu lo one bin. This induces, as noted by IVivekanand et al] 
19981 a periodic long i tude shi ft of successive pulses. Following 
Edwards & Stapners (2003a), we have compensated for this 
longitude shift of each pulse, and thereby avoiding artificial 
features appearing in the spectra that are derived from the 
pulse-stacks. All pulse longitudes in this paper have an arbi- 
trary offset because absolute alignment was not necessary for 
our analysis. 

In the left panel of Fig. [2 one can see a sequence of 100 
pulses of one of the new drifters we have found which clearly 
shows the drifting phenomenon. Drifting means that the sub- 
pulses drift in longitude from pulse to pulse and thereby the 
pulsar emission shows diagonal intensity bands in the pulse- 
stack (drift bands). The drift bands are characterized by two 
numbers: the horizontal separation between them in pulse lon- 

3 http : //pulsar . princeton . edu/tempo/ 



gitude (Pi) and the vertical separation in pulse periods (P3). 
The drift bands of this pulsar are clearly seen by eye in the 
pulse-stack and the values P2 and P3 could in principle be 
measured directly, but in many cases of the newly discovered 
drifters the drift bands are not visible to the eye. To be able 
to detect the drifting phenomenon in as many pulsars as possi- 
ble, all the pulse-stacks were analyzed in a systematic way as 
described in the next two subsections. 

2.3. Processing of the pulse-stacks 

In Fig. [2 the products of our method of analysis are shown for 
three pulsars and in this section it is explained how these plots 
are generated from the pulse-stack and how one can interpret 
them. 

The first thing that is produced from the pulse-stack is the 
integrated pulse profile. This is simply done by vertically inte- 
grating the pulse-stack, i.e. adding the bins with the same pulse 
longitude in the successive pulses: 

1 N ^ 

«=^2>« 0) 

7=0 

Here yu,- is the average intensity at longitude bin i, and S y is the 
signal at pulse longitude bin i and pulse number j in the pulse- 
stack and N is the number of pulses. In Fig. [2 the solid line 
in the top panels corresponds to the integrated pulse profile 
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(which is normalized to the peak intensity). On the horizontal 
axis is the pulse longitude in degrees and the value can be read 
from the horizontal axis of the panel below which is aligned 
with the top panel. 

The first basic method to find out if there is some kind of 
subpulse modulation is to calculate the longitude-resolved vari- 
ance cP: 

1 

j=o 

and the longitude-resolved modulation index m, 

mi = — (5) 
Mi 

The modulation index is a measure of the factor by which the 
intensity varies from pulse to pulse and could therefore be an 
indication for the presence of subpulses. In Fig.^the open cir- 
cles in the top panel is the longitude-resolved standard devi- 
ation <x, and the solid line with error bars corresponds to the 
longitude-resolved modulation index m,. 

The detection of a modulation index does not give infor- 
mation about whether the subpulse modulation occurs in a 
systematic or a disordered fashion. The first step in detect- 
ing a regular intensity variation is to calc ulate the Lon gitude 
Resolved Fluctuation Spectrum (LRFS; Backed Il970n . The 
pulse-stack is divided into blocks of 5 12 successive pulses 4 and 
the Discrete Fourier Transform (DFT) was performed on these 
bloc ks to calculate the LRFS (for details of the analysis we re- 
fer to lEdwards & StapperJ20ci3l2003al) . The fluctuation power 
spectra of the different blocks were then averaged to obtain the 
final spectrum. 

In Fig. ^ the LRFS of the three pulsars are shown below 
the pulse profile plots. The units of the vertical axis are in 
cycles per period (cpp), which corresponds to P0/P3 in the 
case of drifting (where P3 is the vertical drift band separation). 
The horizontal axis is the pulse longitude in degrees, which is 
aligned with the plot above. The power in the LRFS is hori- 
zontally integrated, producing the side panel. If the emission 
of the pulsar is modulated with a period P3, then a distinct re- 
gion of the LRFS will show an excess of power (i.e. a feature) 
in the corresponding pulse longitude range. The LRFS can be 
used to see at which pulse longitudes the pulsar shows sub- 
pulse modulation and with which periodicities. The grayscale 
in the LRFS corresponds to the power spectral density. Under 
Pars eval's theorem, the sum med LRFS is identical to Eq. 4 
( Edwards & Stappers 2003a), so integrating the LRFS verti- 
cally gives the longitude resolved variance (the open dots in 
the plot above the LRFS). 

The detection of a modulation index suggests that there is 
subpulse modulation and by analyzing the LRFS it can be de- 
termined if this modulation is disordered or (quasi-)periodic. 
However from the LRFS one cannot determine if the sub- 
pulses are drifting over a certain longitude range, because 
to calculate the LRFS only DFTs along vertical lines in the 

4 For a few observations with a low number of pulses shorter trans- 
forms were used. 



pulse-stack are performed. To determine if the subpulses are 
drift ing, the Two-Dimensi onal Fluctuation Spectrum (2DFS; 
Edw ards & Stappersl2002l) is calculated. The procedure is sim- 
ilar to calculating the LRFS, but now we select one or more 
pulse longitude ranges between which the DFT is not only cal- 
culated along vertical lines, but along lines with various slopes. 
The effect is that the pulse longitude information that we had 
in the LRFS is lost, but we gain the sensitivity to detect pe- 
riodic subpulse modulation in the horizontal direction (i.e. if 
there also exists a preferred P2 value). Following the same pro- 
cedure used while calculating the LRFS, the pulse-stack is di- 
vided in blocks of 5 12 successive pulses 4 and the spectra of the 
different blocks were then averaged to obtain the final spectra. 

In Fig.[T]the 2DFS is plotted below the LRFS. The vertical 
axis has the same units as the LRFS, but now the units of the 
horizontal axis are also cycles per period, which corresponds 
to P0/P2 in the case of drifting (where P2 is the horizontal 
drift band separation in time units). The power in the 2DFS 
is horizontally and (between the dashed lines) vertically inte- 
grated, producing the side and bottom panels in Fig.^ These 
panels are only produced to make it easier to see by eye what 
the structure of the feature is. 

From the pulse-stack in Fig. [2 one can see that two suc- 
cessive drift bands of PSR B1819-22 are vertically separated 
by P 3 * I8.OP0 and horizontally by P 2 ~ 0.025P - Instead 
of measuring drifting directly from the pulse-stack, we use 
the 2DFS. From both the 2DFS and LRFS of this pulsar we 
see that there are multiple drift features. This is because PSR 
B1819-22 is a drift mode changer (i.e. the drift bands have 
different slopes in different parts of the observation). We note 
that only one drift mode is seen in the short stretch of pulses 
shown in the pulse-stack in Fig.[Q For PSR B1819-22 one can 
see the main feature in the LRFS around 0.056 cpp, which cor- 
responds to the Pj value we see in the plotted pulse-stack. In 
the 2DFS of this pulsar we see the main feature at the same 
vertical position as in the LRFS (corresponding to the same P3 
value) and because the feature is offset from the vertical axis we 
know that the subpulses drift. From the horizontal position of 
the feature in the 2DFS we see that P 2 ~ -Pq/40 = -0.025P , 
which corresponds well with the P 2 measured directly from the 
pulse-stack shown. 

In this paper we use the convention that P3 is always a pos- 
itive number and P 2 can be either positive or negative. A neg- 
ative value of P 2 means that the subpulses appear earlier in 
successive pulses, which is called negative drifting in the lit- 
erature. The tabulated signs of P 2 in this paper therefore cor- 
respond to the drift direction, such that a positive sign corre- 
sponds to positive drifting. To comply with this convention, 
all the plotted 2DFSs in this paper are in fact flipped about 
the vertical axis compared with the definition of the 2DFS in 
lEdwardsfe Starts] ^20021). 

2.3.1 . Interference 

To reduce the effect of interference on the LRFS and 2DFS the 
spectrum of the off-pulse noise was subtracted from the LRFS 
and 2DFS if a large enough off-pulse longitude interval was 
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available. Interference will in general not be perfectly removed 
by this procedure, however any artificial features produced by 
interference can easily be identified because it will not be con- 
fined to a specific pulse longitude range. In Fig. ^ the spectra 
of PSR B2043-04 shows interference with a P /P 3 0.372. 
In the spectra as shown in appendix[X]and||}] the channels con- 
taining interference are set to zero, thereby improving the vi- 
sual contrast of the plots. 

In this paper the modulation index is not directly derived 
from the pulse-stack (Eqs.|4]and|5}, but from the LRFS. This is 
done by vertically integrating the LRFS, which gives the lon- 
gitude resolved variance (Eq.|4j. The advantage of this method 
is that by excluding the lowest frequency bin the effect of inter- 
stellar scintillation (which at this observing frequency has typ- 
ical low frequencies) can be removed from the modulati on in- 
dex (for deta ils of the analysis we refer to lEdwards & S tap pers 
l2002tl2003all . After exclusion of the lowest frequency bins the 
variance is overestimated by m 2 . ,. where m 2 . , is the mod- 

J scint' scmt 

ulati on index induced by the scintillation (see Eqs. 20-22 of 
lEdwards & Stappersl2004l) . The longitude resolved modulation 
index and standard deviation are corrected accordingly. 

2.4. Analysis of the drift features 

If a feature is seen in the 2DFS and we make sure it is not asso- 
ciated with interference, Pi and Pj can be measured and its sig- 
nificance determined. The drift feature will always be smeared 
out over a region in the 2DFS. This could be because there is 
not one fixed value of the drift rate throughout the observation 
due to random slope variations of the drift bands, drift mode 
changes or nulling. But the feature is also broadened if the drift- 
ing is not linear (i.e. subpulse phase step s or swings) and be- 
cause of subpulse amplitude windowing ( Edw ards & Stap pers 
120021) . 

Because of all these effects it is impossible to fit one spe- 
cific mathematical function to all the detected features, so it 
is more practical to calculate the centroid of a rectangular re- 
gion in the 2DFS containing the feature. The advantage of this 
procedure is that no particular shape of the feature has to be 
assumed. The centroids are defined as: 

Po/Pi = jYj y,Fu 

k,l 

P0/P2 = jYj xkFu (6) 

k,l 

Here k and I are indices for the horizontal and vertical position 
within the region in the 2DFS containing the feature, xt, and yi 
are the corresponding axis values of bin (k, I), Fti is the power 
in that bin and F = Y*k ,1 Fid is the total power in the selected re- 
gion of the 2DFS. An uncertainty on the position of the centroid 
(hence on Pi and P3) can be estimated by using the power in a 
region in the 2DFS that only contains noise. However we have 
found that in most cases the uncertainty is dominated by the 
more or less arbitrary choice of what exact region in the 2DFS 
is selected around the feature. To estimate this uncertainty we 
have therefore calculated the centroid for slightly different re- 
gions containing the feature. 



Because the analyzed pulse sequences are sometimes rel- 
atively short, there is the possibility that the occasional oc- 
currence of strong (sub)pulses are dominating the spectra and 
therefore lead to misleading conclusions. To estimate what kind 
of "random" fluctuation features one can expect from a given 
pulse sequence, we have randomized the order of the pulses and 
then passed this new pulse-stack through our software to deter- 
mine the magnitude of any features which could be attributed 
to strong (sub)pulses. Any actual drifting will lose coherence in 
this process and thus we can use the randomized results when 
estimating the significance of drift features when there is no 
well defined P3 . 

For most pulsars there is power along the vertical axis in the 
2DFS. Therefore the centroid of a larger region around a drift 
feature usually results in a centroid located closer to the vertical 
axis, hence in a larger absolute value of Pi. In most cases this 
causes the uncertainty on Pi to be asymmetric around its most 
likely value and therefore it is useful to tabulate the uncertainty 
in both signs. 

If the centroid of the feature in the 2DFS is significantly 
offset from the vertical axis, it means that Pol Pi + and that 
drifting can be associated with the feature (i.e. there exists a 
preferred drift direction of the subpulses from pulse to pulse). 
Drifters are defined in this paper as pulsars which have at least 
one feature in its 2DFS which show a significant finite Pi. It 
must be noted that drift bands are often very non-linear and 
therefore the magnitude of Pi is probably of little meaning. 
Also we want to emphasize that not only an offset from the 
vertical axis indicates drifting, but any asymmetry about the 
vertical axis indicates drifting-like behavior. Any frequency de- 
pendence along the vertical axis indicates structured subpulse 
modulation, either quasi-periodic or as a low frequency excess. 
Pulsars with a low frequency excess show an excess of power 
in their spectra toward long frequencies (i.e. the spectra are 
"red"). 

The drifters are classified into three classes and an exam- 
ple of each class is shown in Fig. One can see that PSR 
B2043-04 has a vertically narrow drift feature in its spectra, 
meaning that Pj has a stable and fixed value throughout the 
observation. We will call these pulsars the coherent drifters 
class (class Coh in table 0. The criterion used for this class 
is that the drift feature has a vertical extension smaller than 
0.05 cpp. The pulsars that show a vertical, broad, diffuse drift 
feature are divided into two classes, depending on whether the 
feature is clearly separated from the alias borders (P0/P3 = 
and P0/P3 = 0.5). In table[2]the pulsars in class Dif are clearly 
separated from the alias border and the pulsars in class Dif* 
not. In Fig.Q]PSR B1819-22 is a diffuse (Dif) drifter and PSR 
B2148+63 is a Dif drifter. 

Besides the drifters there is also a class of pulsars which 
show longitude stationary subpulse modulation (class Lon in 
table |3. These pulsars show subpulse modulation with a P3 
value, but without a finite Pi value. Because it is not clear if 
these pulsars should be counted as drifters or as non-drifters, 
they are excluded from the statistics. 

For many pulsars we find that the magnitude of Pi ex- 
ceeds the pulse width. This means in the case of regular drift- 
ing that in a single pulse only one subpulse is visible and that 




Fig. 2. A section of the pulse-stacks and the derived spectra of the two artificially generated pulse sequences to illustrate possible 
effects of pulsars that constantly change their alias order. The drifting in the left pulse-stack constantly crosses the P3 = 2Pq alias 
border and the right pulse-stack constantly changes its apparent drift direction via longitude stationary subpulse modulation. For 
the explanation of the plots we refer to Fig.[Jand the main text. 



the drifting will manifest itself more as an amplitude modu- 
lation rather than as a phase modulation. An illustrative ex- 
ampl e of a reg ular drifter with a large P2 is PSR B0834+06 
( Asgek ar & DeshpandeH2005l) . Whether the drifting manifests 
itself more as an amplitude modulation or as a phase modu- 
lation will largely depend on the viewing geometry. Also the 
presence of pulse sequences without organized drifting, longi- 
tude stationary subpulse modulation or drift reversals will re- 
sult in a large P% value. 

It must be noted that the calculation of the 2DFS is an av- 
eraging process. This is what makes it a powerful tool to de- 
tect drifting subulses in low S/N observations, but at the same 
time this implies that different pulse-stacks can produce similar 
2DFS. For instance a feature that is split and shows a horizontal 
separation can be caused by drift reversals, but al so by subpulse 
phase steps or swings (see Edw ards et alj |2003 for a pulse se- 
quence of PSR B0320+39 and the resulting 2DFS). Note also 
that with only the LRFS it is impossible to identify complex 
drift behavior like drift reversals or subpulse phase steps. 

The details of each observation can be found in table |2] 
including the classification we made, the measured P2 and P3 
values, the modulation index and the detection threshold for 
the modulation index. It must be emphasized that the P2 and 
P3 values are average values during the observation. Especially 
when the pulsar is a drift mode changer a different observation 
may lead to different values for Po and P3. The plots of all the 



pulsars can be found in appendix [X] For some of the pulsars 
the 2DFS for two different pulse longitude ranges are shown 
if useful. The plots of these pulsars come after the plots of the 
pulsars with only one 2DFS plot. The same plots can also be 
found in appendix|B] but there they are ordered by appearance 
in the text instead of ordered by name. 

2.4.1 . Drift reversals 

As can be seen from Fig. [2 the subpulses of PSR B1819-22 
appear to drift toward the lea ding part of the pulse profile. I n 
the sparking gap model (e.g. iRuderman & Sutherlandin"975h . 
every subpulse is associated with one emission entity close to 
the surface of the star. These entities (the sparks) move around 
the magnetic axis, causing the subpulses to drift through the 
pulse window. Because the emission entities are only sampled 
once per rotation period of the star, it is very difficult to de- 
termine if the subpulses in one drift band correspond to the 
same emission entity for successive pulses. For instance for 
PSR B18 19-22 we do not know if the emission entities drift 
slowly toward the leading part of the pulse profile (not aliased) 
or faster toward the trailing part of the pulse profile (aliased). 

The physical conditions of the pulsar probably determines 
what the physical drift rate of the emission entities are, rather 
than the observed (possibly aliased) drift rate. This could al- 
ready be a serious problem if one wants to correlate physi- 
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cal parameters of pulsars to the observed drift rate of coher- 
ent drifters, for which we at least know they stay in the same 
alias mode. If a feature in the 2DFS is not clearly separated 
from the alias borders, the power in that feature could consist 
of drift bands in different alias modes (i.e. the apparent drift di- 
rection could be changing constantly during the observation). 
In that case the measured value of P3 using the centroid of the 
feature is related to the true drift rate of the emission entities 
in a complicated way, depending on what fraction of the time 
the pulsar spends in which alias mode. Therefore it is expected 
that it will be very hard to find a correlation between physical 
pulsar parameters and the P3 values of the pulsars in the Dif* 
class, so it will be useful to classify the pulsars depending on 
whether the features in the 2DFS are clearly separated from the 
alias borders. Inspection of the pulse-stacks with strong enough 
single pulses reveals that some of the pulsars in the Dif* class 
change their drift direction during the observation. This is fur- 
ther evidence to indicate the value of considering the Dif and 
Dif* classes separately. 

To illustrate this we have artificially generated two pulse 
sequences of a pulsar that has a variable rotation period of the 
emission entities in two different scenarios (see Fig.|5Jl. In both 
scenarios the emission entities are simulated to drift from the 
trailing to the leading edge of the pulse profile with a variable 
drift rate. In the left sequence the vertical drift band separa- 
tion P3 is close to 2fo and in the right sequence the P3 period 
is much larger. In the left sequence the subpulses around the 
first pulse appear to drift toward the leading edge of the pulse 
profile. As time increases, we speed up the rotation of the emis- 
sion entities, which causes P3 to become smaller. Around pulse 
number 15 the Nyquist border P3 = 2Pq is reached and the drift 
pattern becomes a check-board like pattern. As time further in- 
creases the emission entities are still speeding up, causing clear 
drift bands to reappear with an opposite apparent drift direc- 
tion (around pulse 25). After this the rotation of the emission 
entities is gradually slowed down to the initial value, causing 
the drift bands to change apparent drift direction again around 
pulse 50. The same cycle is repeated for the next pulses. In this 
simulation the carousel rotation period is set to vary with about 
40% around its mean value. 

The resulting spectra of this pulse sequence are also shown 
in Fig. 13 The LRFS shows that the subpulse modulation is ex- 
tended toward the P3 = 2Pq alias border and the 2DFS shows 
a feature that is split by the vertical axis, because there are 
two apparent drift directions in the pulse sequence. As can be 
seen in the bottom panel, there is more power in the left peak. 
This corresponds to more power being associated with negative 
drifting (drifting toward the leading edge of the pulse profile). 
This is also directly visible in the pulse-stack. A good exam- 
ple of a known pulsar t hat shows thi s kind of drift behavior is 
PS R B230 3+30 fe.g. lRedman et alJ E"005) and its 2DFS (see 
Fig. IA.15l indeed shows a very clear double peaked feature. 

In the second scenario of Fig.|2]the pulse-stack shows drift 
bands with a much larger Pj. Because of possible aliasing this 
does not directly imply that the emission entities are rotating 
slower. In fact, we have chosen the entities to rotate faster than 
in the first scenario. This causes the driftbands to be aliased and 
the drift bands appear to drift in the opposite direction to the 



PSR B0809+74 PSR B0809+74 




Fig. 3. This figure shows the spectra of the well known drifter 
PSR B0809+74 (left) as well as the results from the same data 
after putting the pulses in a random order (right). 



emission entities (which again are simulated to drift toward the 
leading edge of the pulse profile). As time increases the rota- 
tion of the emission entities is sped up. Because of aliasing the 
drift rate appears to decrease until around pulse 25 the drifting 
has become longitude stationary (P3 — 00). The emission enti- 
ties are now rotating so fast that in one pulse period time they 
exactly reappear at the pulse longitude of another drift band. 
When the rotation period of the emission entities is speeded up 
even further, the drift bands are changing their apparent drift 
direction again as can be seen around pulse 50. Now the rota- 
tion period of the emission entities is slowed down again until 
around pulse 100 the initial conditions are reached again. After 
this the same cycle is repeated. A clear example of a known 
pulsar tha t shows this kind of drift reve rsals is PSR B0826-34 
(Gunta et alfcnofllEsamdin et al J200.4. 



The LRFS of this sequence shows that the subpulse mod- 
ulation is extended toward the horizontal axis and the 2DFS 
shows again a feature that is split by the vertical axis. As can 
be seen in the bottom panel, there is more power in the right 
peak, which corresponds to more power associated with pos- 
itive drifting. In the pulse-stack there are indeed more drift 
bands that drift toward the trailing edge of the pulse profile 
than in the opposite direction. 
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2.4.2. Non linear drift bands 

In the left panel of Fig. [5] the spectra of the well known reg- 
ular drifter PSR B0809+74 are shown. The drift feature in 
the 2DFS shows a clear horizontal structure caused by the 
non linear dri ft bands of PSR B0809+ 74 at this observing fre- 
quen cy (Edwards & Stappers 2003b; Proszvriski & Wolszczanl 
1986; Wols zczan et alJfl98 lh . If drift bands are not straight, 
there is no one unique value of P2 that is associated with the 
drifting and therefore the drift feature in the 2DFS will be more 
complex than just one peak. This makes P2 an ill-defined pa- 
rameter. However this is no shortcoming of the 2DFS, it is a 
shortcoming of the whole concept of P2 under curved drift- 
bands. In this paper P2 is just a rough measure of the presence 
of drift, its direction and the magnitude of the slope in a overall 
mean sense only. 

If the drift bands are non linear, the subpulses will have a 
pulse longitude dependent spacing. This pulse longitude depen- 
dent spacing is described by the so-called "modulation phase 
profile" or "phase envelope" fedwards & Stappersl2002l) . 

2.4.3. Quasiperiodic subpulses 

For most drifters the vertical drift band separation P3 is a much 
better defined parameter than the horizontal separation P2, e.g. 
the drift feature in the 2DFS of PSR B0809+74 (left panel of 
Fig-El is much sharper in the vertical direction than in the hor- 
izontal direction. There is however the possibility that single 
pulses show regular spaced subpulses, while there is no mem- 
ory for where the subpulses appear in successive pulses. In that 
case there is a P2 value, but P3 is undefined. One can simu- 
late such a scenario by putting the pulses of a regular drifter in 
a random order. This is done for PSR B0809+74 in the right 
panel of Fig.[3] 

One can see, first of all, that the longitude-resolved standard 
deviation and modulation index are independent on the order- 
ing of the pulses, as expected from Eqs.|4]and[5] Secondly, the 
spectra do not show any vertical structure anymore. This indi- 
cates that, as expected, P3 has become undefined. The 2DFS is 
symmetric about the vertical axis, so there is no preferred drift 
direction anymore. The horizontal subpulse separation is how- 
ever still visible in the 2DFS as two vertical bands at +37 cpp, 
the same horizontal position as the largest peak in the 2DFS in 
the left panel of Fig. [5] 

We have found two pulsars that shows these kind of fea- 
tures: PSRs B2217+47 and B0144+59. In the carousel model 
this phenomenon could be explained by a highly variable cir- 
culation time that causes the alias order to change constantly. 
However there is no evidence that this phenomenon is related 
to the same origin as the drifting subpulses, as it only shows 
that the subpulses appear quasiperiodic. 

3. Individual detections 

3.1. Coherent drifters (class "Coh") 

The coherent drifters are the pulsars which show a drift fea- 
ture in their 2DFS over a small P3 range (smaller than 0.05 



cpp). First the pulsars that were already known to have drifting 
subpulses are described followed by the pulsars that were not 
known to show this phenomenon. 

3.1 .1 . Known drifters 

B0148-06: Both components of the pulse profile of this 
pulsar 5 have the same drift sign and P3 varies slightly during 
the observation (see Fig. \A.1H . The drift bands are clearly 
visible in the pulse-stack and the drifting is most prominent in 
the leading component. The feature in the leading component 
seems to consist of different drift modes. Th is is all consistent 
with results reported by IBiggs et al] dl985h . who discovered 
the drifting subpulses at 645 MHz. 

B0320+39: This pulsar is known to show very regular drifting 
subpulses dlzvekova et all [l982), which is confirmed by 
the very narrow drif t feature in our observation (Fig. IA. II . 
Ilzvekova et all ll 1993b have shown that drifting at both 102 and 
406 MHz occurs in two distinct pulse longitude intervals and 
that the energy contribution in the drifting subpulses is less at 
higher frequencies, especially in the leading component. We 
see that this trend continues at higher frequencies, because the 
feature is most prominent in the trailing component at 1380 
MHz. The drift feature in the 2DFS shows a clear horizontal 
structure caused by a subpulse phase step in the drift bands (as 
is also seen for instance for PSR B0809+74 at this frequency). 
At both 328 MHz and 1380 MHz the drift bands of PSR 
B032 0+3 9 are known to show a ph ase step fedwards et alJ 
2003 and Edwards & Stappers 2003b respectively). This is the 
same observation as used by[]Edwards & S tappers! J2003bl) . 
B0809+74: The drift feature in the 2DFS (Fig. |a3J» shows 
clearly horizontal structure, which is caused by a phase step of 
the drifting subpulses. This phase step is only present at high 
frequenci es and is consistent with IWolszczan et alJ dl98ll) at 
1.7 Ghz. IProszvfiski & Wolszczanl dl986l) at 1420 MHz and 
lEdwards & Stappers! J2003bl) at 1380 MHz f this is the same 
data as used by Edwards & S^apper^|2QQ3b|). The drift rate 
is affected by nulls dTavlor & Hugueninlll97ll) and deta iled 
analys is of this phenomenon allowed Ivan Leeuwen et all 
(2003) to conclude that the drift is not aliased for this pulsar. 
B0818-13: This pulsar has a clear drift feature (Fig.lTOt that 
contains almost all power in the 2DFS. The drift feature has a 
horizontal structure like observed for PSR B0809+74 and PSR 
B0320+39. From the modulation phase profile it follows that 
the drift bands make a smooth swing of about 200° in subpulse 
phase in the middle of the pulse profile. A decrease of the drift 
rate in th e middle of the pu lse profile has been reported at 645 
MHz by Big gs et all (U287) and this effect appears to be more 
pronounced in our observation at 1380 MHz. The longitude 
resolved modulation index shows a minimum at the pos ition 
of th e subpulse phase swing (consistent with IBi ggs et all 
l!987h . something that is also observed for the phase steps of 
PSR B0320+39 and B0809+74. The subpulse phase swing 
is clearly visible in the pulse-stack of this pulsar. This phase 

This pulsar is not in our source list, because the required obser- 
vation length is too long. Therefore this pulsar is not included in the 
statistics. 
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swing seems to be related to the complex polarization behavior 
of the single pulses as observed by Edwards (2004). For this 
pulsar nullin g wa s show n to interact with the drift rate by 
ILvne & AshworthI (f l983) and from d e tailed analysis of this 
phenomenon Ijanssen & van LeeuwenI d2004l) concluded that 
the drift of this pulsar is aliased. 

B1540-06: The 2DFS of the two halves of the pulse profile 
show opposite drift directions (Fig. IA.20> . The drift sense 
of the leading part of the pulse profile is consistent with 
the positive drif ting of this pulsar as has been reported by 
lAshworfhl i 19821) at 400 MHz. 

B2045-16: Only the 2DFS of the trailing component is plotted 
in Fig. IA. 141 because we do not detect any features in the 
leading component. However, dri fting has been obse rved 
previously in both components (e.g. lOster & Sieberll 19771) and 
is observed to b e broad with P3 va lues between 2 and 3Pp 
llOster & Sieberlll977l at 17 20 MHz. iNowakowski et"aTH "l982 
at at 1.4 and 2.7 GHz and iTavlor & Hugueninlll97ll at low 
frequencies). The second component in our observation has a 
clear narrow drift feature. This could indicate that this pulsar 
is a drift mode changer and that our observation was too short 6 
to detect the drift rate variations. 

B2303+30: This pulsar is know n to drift close to the alias 
border (e.g. lSieber & Osterll975l at 430 MHz). This pulsar has 
a clear double-peaked feature in its 2DFS exactly at the alias 
border (Fig. IA.15l . which suggest that the apparent drift sense 
changes during the observation because the alias changes 
during the observation. The change of drift sense can be seen 
b y eye in the pulse-sta ck and also in the pulse-stacks shown 
in Redm anet alJ J2005). These authors show that, besides this 
Po,x2Pq 'B' drift mode, there is occasionally also a Pj,~3Po 
'2' drift mode if the S/N conditions are good. There is no 
evidence for this weaker 'Q'-mode in our observation. 
B2310+42: The two components of this pulsar are clearly 
drifting at the alias border (Fig. IA.23> . The drift feature 
in the trailing component is clearly double peaked, so the 
alias mode is probably changing during the observation. 
The leading component shows the same feature at the alias 
border, although much weaker. The domin ant drift sense is 
consistent with the positive drifting found b v I Ash worthl (H 9 8 2 ) 
at 400 MHz. We also find that the low frequency excess of 
the leading component is clearly drifting with two signs as 
well (P 2 = 65° + 35, P 3 = 21 ± 3P and P 2 = -65° ± 15, 
P3 = 70 + lOPoX but no significant drift is detected in the low 
frequency excess of the trailing component. 
B2319+60: This pulsar shows a clear drift component in the 
2DFS of the trailing component of the pulse profile and a 
less clear drift component in the leading component (Fig. 
IA.23> . No significant drifting is detected in the 2DFS of 
the center part of the pulse profile (which is not plotted). It 
has been found that this pulsar is a drift mode changer and 
that the allowed drift m ode transitions follow certain rules 
JWri ght&Fowlerlll98ll at 1415 MHz). In the 2DFS there is 
only evidence for one stable drift mode, but drifting is detected 



6 Our observation was too short to contain the required minimum 
of one thousand pulses. Because we have detected drifting, we did not 
reobserve this pulsar. 



over the whole P3 range. Drift mode changes can been seen 
in the pulse-stack. The drift mode w e see in the 2DFS is the 
'A' drift mode, which was found bv lWright & Fowlerl J 1 98 lit 
to be the most common drift mode. The many nulls in the 
pulse-stack probably smears out the feature over a large P3 
range. There is a sharp P3 130Po feature in the second 
component, which could be related to nulling or mode changes. 



3.1 .2. New drifters (class "Coh") 

B0149-16: The 2DFS of this pulsar (Fig. |A7TJ shows a weak 
drift feature. 

B0609+37: Almost all power in the 2DFS of this pulsar is in 
a well confined drifting feature (Fig. IA.2L meaning that the 
subpulse drifting is very organized and stable. 
B0621-04: A strong and very coherent P3 =s 2Pq feature is 
seen in the LRFS and 2DFS of this pulsar (Fig. IA.2> . which 
is also seen in other observations. The feature is significantly 
offset from the vertical axis, so this pulsar shows very stable 
drifting subpulses. Only the 2DFS of the trailing peak is 
plotted. 

J1650-1654: This pulsar shows a very clear drift feature in 
its 2DFS (Fig. \A.5i . The feature seems to show some vertical 
structure, which could be because of drift rate variations. The 
feature also seems to show some horizontal structure like 
PSR B0818-13, which could indicate that the drift bands are 
curved or show a phase step. This would be consistent with the 
minimum in the modulation index in the middle of the pulse 
profile. However this observation is too short to state if this 
effect is due to a systematic drift rate change across the pulse 
profile or due to random variations. 

B1702-19: The pu lse profile of this pulsar shows an interpulse 
<Biggs_elalj[l988). The main pulse shows a drift feature 
and the 2DFS of the interpulse shows a feature with the 
same P3 value (Fig. IA.20> . but no significant offset from the 
vertical axis could be detected for this feature. This is not 
the only pulsar to show correlations in emission properties 
between the main- and interpulse. PSR B 1822-09 exhibits 
an anti-correlation bet ween the intensity in the main- and the 
interpulse (e.g. lFowler et al.ll98ltlFowler & Wrighl l982) and 
also a correlation in the subpulse modulation (e.g. iGil et alJ 
Il994l) . That pulsar also shows drifting in the main pulse in 
the 'fi'-mode and longitude stationary subpulse modulation in 
both the main- and interpulse in the '2' -mode with the same 
P3. Also PSR B 1055-52 is known to show a main pulse- 
interpulse correlation. For that pulsar a correlation between the 
intensi ty of the main and interpulse has been found by Biggs 
dl990l) . However PSR B 1702- 19 is the first pulsar that shows 
a correlation between the drifting subpulses in the main pulse 
and the subpulse modulation in the interpulse. 
B1717-29: A very narrow drift feature is seen in both the 
LRFS and the 2DFS of this pulsar (Fig. IA.5t . Because the very 
low S/N of this observation it was impossible to measure a 
significant modulation index using the whole P3 range of the 
LRFS. By only using the frequencies in the LRFS that contains 
the drift feature it was possible to significantly determine 
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the modulation index corresponding to the drift feature. This 
drifting is confirmed in another observation. 
B1839-04: Both components of the pulse profile of this 
pulsar are drifting (Fig. IA.2U . The drift bands are clearly 
visible in the pulse-stack to the eye and both components have 
an opposite drift sense. The slope of the drift bands in the 
two components are mirrored and the drift bands of the two 
components are also roughly in phase. So when a drift band is 
visible in one component, it is also visible (although mirrored) 
in the other component. This "bi-drifting" subpulse behavior is 
also observed for PSR J08 15+09, which also has opposit e drift 
senses in different components JlVIcLaughlin et alJ |2004). This 
"bi-drifting" could be a sign that these puls ars have bot h an 
inner annular gap and an inner core gap rtOiao et all2004l). but 
also "double imaging" could be responsible ( Edwards et alJ 
120031) . Note also that the second harmonic of the drift feature 
is visible in the 2DFS of especially the trailing component. 
B1841-04: This pulsar has a weak, definite drift feature in its 
2DFS (Fig.lA!9t. which is also visible in the LRFS. 
B1844-04: There is a weak detection of a narrow drift feature 
in the 2DFS of this pulsar, which is also visible in the LRFS 

(Fig.EIB- 

J1901-0906: The trailing component of this pulsar shows 
a clear and narrow drift feature in its 2DFS, which is not 
detected in the leading component (Fig. IA.2H . The 2DFS of 
the leading component has a drift feature with a different P 3 
which is also present in the right component (P2 = -8.2° ± 9.2, 
P3 = 6.8 + 0.5Po). The drifting can be seen by eye in the 
pulse-stack. The different measured P3 values in the two 
components could indicate that this pulsar is a drift mode 
changer. 

B2000+40: Although this observation is contaminated by 
interference, clear drifting is detected in the leading com- 
ponent. The rest of the pulse profile (mostly in the trailing 
component) is also drifting (Fig. IA.22> . The feature in the 
leading component shows horizontal structure which could be 
caused by drift reversals or more likely by a subpulse phase 
jump or swing. 

B2043-04: This pulsar has a very clear and narrow drifting 
component in its 2DFS (Fig. IA. 141 . The feature is perhaps 
extended toward the alias border, but this is not significant. 
Almost all power in the 2DFS is in the drift feature. 



3.2. Diffuse drifters (classes "Dif" and "Dif* ") 

These pulsars show a drift feature over an extended Pj range 
(larger than 0.05 cpp). If the drift feature is clearly separated 
from both alias borders (P Q /Pi = and P Q /P 3 = 0.5), the 
pulsar is classified as Dif. However if it is not the pulsar is 
classified as Dif* in table |2 In this section the pulsars in the 
latter class are indicated with an asterisk next to their name. 
Note that not all drift features in the spectra have peaks which 
are offset from the vertical axis, but they must be asymmetric 
about the vertical axis. 



3.2.1 . Known drifters 

B0031-07: This pulsar shows a broad drifting feature in its 
2DFS (Fig. IA.lt . Three d rift mo des have been found for this 
pulsar by Hugueni n et alJ at 145 and 400 MHz. In our 

observation most power in the 2DFS is due to the 'A' -mode 
drift (P3 = 12Po). The slope of the drift bands change from 
band to band (e.g. IVivekanand & JoshU Il997h . causing the 
feature to extend vertically in the 2DFS. The 'B'-mode drift 
(P3 = 6P0) is also visible in our observation, but there is no 
feature corresponding to 'C'-mode drift ( P3 = 4Pp). This i s 
consistent with the multifrequency study of Smit s et al .(2005). 
B0301+19*: The trailing component shows a broad drift 
feature in its 2DFS (Fig. IA. 171 . but no drifting is detected 
in the leading component. This pulsar is observed to have 
straight drift bands in both components of the pulse-stack 
dSchonhardt & Sieberlfl973l at 430 MHz). The feature in the 
trailing comp onent is reported to be broader than in the leading 
component (backer et al.lll975l also at 430 MHz), probably 
because drifting subpulses appear more erratic in the trailing 
component. The feature we see is also broad and may even be 
extended to the alias border. 

B0329+54*: The power in the LRFS of this pulsar peaks 
toward I/P3 = 0, as reported by iTavlor & Hugueninl dl 97 ll) 
for low frequencies (Fig. IA.17l . Drifting is detected in four of 
the five components. The third component (the right part of 
the central peak) has a broad drift feature and the subpulses 
have a posit i ve dri ft sense, something that is also reported by 
iTavloret aDd 19751) at 400 MHz. Besides these known features 
we find that the first component (left peak) and the fourth 
component (the bump between the central and trailing peak) 
are also drifting with a positive drift senses: P2 = 750° ± 900 
and 147° ± 1 10 and P 3 = 2.9 + 1.7 and 4.2 ± 1.0 respectively. 
The second component (the left part of the central peak) has 
an opposite drift sense: P 2 = -175° + 150 and P 3 = 3.0 ± 1.7. 
The last component shows no significant drifting. The 2DFS 
of the second and third component are shown in Fig. IA. 171 
The difference between the values of P3 in the different 
components seems not to be significant. 
B0628-28*: Sporadic drifting w ith a positive d rift sense has 
been reported for this pulsar bv lAshworthl(ll982l) at 400 MHz, 
but the P2 and P3 values could not be measured. The positive 
drift sense is confirmed in our observation as a clear excess 
of power in the right half of the 2DFS (Fig. IA.3l . The feature 
in the 2DFS is not separated from either alias borders. Most 
power in both the LRFS and 2DFS is in the lower half. 
B0751+32*: The 2DFS of the leading component of the pulse 
profile of this pulsar 5 shows drifting over the whole P3 range 
with a negative Pi value (Fig. IA.19l . This can clearly be seen 
in the bottom plot, which shows an excess of powe r in the 
left ha lf. This confirms the drifting as reported by iBackusl 
(1981) at 430 MHz. Both components also show a strong 
P3 = 70 ± lOPo feature. This feature shows negative drifting 
in the leading component, but no significant drifting in the 
trailing component. 

B0823+26*: Only the pulse longitude range of the main pulse 
is shown in Fig. I A. 3 1 an d the 2DFS of the main pulse shows a 
clear broad drift feature. Backer ( 1973) found that at 606 MHz 
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this pulsar shows drifting in bursts, but the drift direction is 
different for different bursts. In our observation there seems to 
exist a clear preferred subpulse drift direction, so this pulsar is 
classified as a drifter. 

B0834+06*: The 2DFS of both components have a weak 
drift feature at the al ias bord er (Fig. lA.l9| i. This confirms the 
drifting detected by ISutton et alJ dl97Qh at low frequencies. 
The circulation time of this p ulsar (P3) has been measured by 
lAsgekar & Deshpanddll2005l) . 

B1133+16*: The 2DFS of both components of this pulsar 
(Fig. IA.19I I show a very broad drifting feature with the same 
drift sign consistent with other data we have analyzed. The 
trailing component shows also a long period drift feature 
(P 2 = 160° + 100 and P 3 = 33 ± 3P n ). This drifting is 
consistent with the drifting found by Nowakowski ( 1996) at 
430 and 1418 MHz and bv lBackeil ( fl973l) : lTavlor et alJ Jl975h 
at low frequencies. 

B1237+25*: The 2DFS of the outer components of the pulse 
profile are clearly drifting with opposite drift sign (they are 
plotted in Fig. IA.20l i. The 2DFS of the three inner components 
(which are not plotted) all show drifting with a positive drift 
sense (except the middle one which does not show significant 
drifting). The values are P 2 = 26° ± 18, P 3 = 2.7 ± IP and 
P 2 = 17° + 9, P3 = 2.7 + IPn respectively. Th is drifting is 
consistent with IProszvfiski & Wolszczanl Jl986l) at 408 and 
1420 MHz. 

J1518+4904: This millisecond pulsar 7 has a clear broad drift 
feature (Fig. IA.4> . This pulsar wa s already known to be a 
drifter (Edwards & S tappers! Eo03al at 1390 MHz), showing 
that drifting is not an phenomenon exclusive to slow pulsars. 
B1642-03*: Drifting is observ ed to occur in bursts in th is 
pulsar with both d rift senses jTavlor & Hugueninl Il97ll at 
400 MHz) and also iTavlor et alJll975l report that there is no 
preferred drift sense at 400 MHz. The 2DFS of our observation 
(Fig. IA.4> reveals a broad drift feature with a preferred drift 
sense, so this pulsar is classified as a drifter. The alias border 
seems to be crossed on both sides, because the feature is 
extended over the whole P3 range and seems double peaked. 
B1822-09*: For this pulsar a correlation in the subpulse 
modulation between the main- and interpulse has been found 
(see the text of the coherent drifter PSR B 1702- 19 for details). 
There are no features in the spectra of our observation of the 
the interpulse and therefore the interpulse is not plotted in 
Fig. IA.2f1 There is drifting detected in the trailing component 
of the main pulse, but it is not clear what exact range in the 
2DFS shows drifting causing the large uncertainty on the P3 
value. T he observation is co nsistent with the 'fi'-mode drift 
found bv lFowler et ail dl98ll) at 1620 MHz with a P 3 * 1 1P - 
There is also a feature at 0.02 cpp, which cou ld be rel ated 
to the Pj ^ 40Po '<2' -mode drift found by iFowler et alJ 
1J). Contrary to their results, in our observation there is no 
evidence that this feature is drifting. This could be because our 
observation was much shorter. 

B1845-01*: The 2DFS of this pulsar (Fig. lA~TUt shows a 
broad drifting feature confirming the detection of drifting in 



7 This pulsar is not in our source list, because the flux is not in the 
catalog. Therefore this pulsar is not included in the statistics. 



this pulsar bv lHankins & Wolszczanl i 19871) at 1414 MHz. 
B1919+21: Both components of this pulsar are clearly drifting 
and almost all power in the 2DFS is in the drift feature (Fig. 
IA.22I I. The feature of the leading component shows horizontal 
structure. The centroid of the whole feature in the leading 
component gives P 2 = -13° ± 3 with the same P3 value. 
The reason for this horizontal structure in the drift feature is, 
like PSR B0809+74, that there is a subpulse phase step in 
the drift band s. This observation confi rms the reported phase 
step by IProszvfiski & Wolszczanl Jl986h seen at 1420 MHz. 
B1929+10*: The LRFS peaks at low frequencies (Fig.lA~T3l . 
comparable to what wa s found bvlNowakowski et alJ i 19821) at 
0.43, 1.7 and 2.7 GHz. lOsteret alJll97 7 suggested, using 430 
MHz data, that this pulsar drifts. The 2DFS of our observation 
shows two broad features with opposite drift sense with two 
different P3 values. The most clear drift feature is between the 
dashed lines and the oth er feature is dir ectly above this feature 
up to I/P3 ^ 0.3. Also Backer ( 19731) has seen two features 
in the LRFS of this pulsar at 606 MHz and the short period 
feature appeared to have a negative drift and the long period 
fluctuations appeared to be longitude stationary. A negative 
drift sense is detected for the short period feature, but the long 
period feature shows a positive drift. Both drift features are 
arising from the leading half of the pulse profile. An explana- 
tion for the observed behavior is that this pulsar is a drift mode 
changer showing different P3 values with opposite drift senses. 
There is also an indication for a P3 ^ 2Pq modulation. There 
is a strong very narrow spike around P3 = 5Po, which could be 
caused by a few strong pulses. 

B1933+16*: This pulsar shows subpulse modulation o yer the 
whole P3 range (Fig. IA. 13i . It was found by iBackerl dl973h 
that there is no preferred subpulse drift sense at 430 MHz, 
however regular d rifting with P3 ^ 2.2fo has been reported by 
Osteretal I II 19771) at 430 MHz. We can confirm that there is 
preferred positive drifting in a broad feature near the P3 = 2 
alias border. 

B1944+17*: This pulsar shows a clear broad drift feature in 
the 2DFS (Fig. IA. 131 and the drifting can clearly be seen 
by eye in the pulse-stack. The fea ture is broad bec ause this 
pulsar shows drift mode changes jDeich e t al. 1986 at both 
430 and 1420 MHz). The P 3 = 13P 'A'-mode drift and the 
P 3 = 6AP Q 'fl'-mode drift are visible in the 2DFS at 0.075 
and 0.16 cpp respectively. We see also evidence for a feature 
in a different alias mode, although much weaker than the main 
feature (P 2 = 12.4° + 0.7 = 30 cpp an d P3 = 22.2 ± 14P ). It 
could be that the zero drift 'C'-mode JDeich et alj fl986) is a 
drift mode for which the drift sense is changing continuously. 
B2016+28*: The 2DFS of the trailing part of the pulse profile 
shows a very broad drifting feature (Fig. IA.22L which is 
caused by drift mode changes (e.g. lOster et alJll977l at both 
430 and 1720 MHz). The leading part of the pulse profile 
shows the same broad drift feature, but also a much stronger 
slow drift mode. This slow drift mode is probably also seen in 
the trailing part of the pulse profile, but less pronounced. The 
drift bands can be seen by eye in the pulse-stack. 
B2020+28*: The LRFS shows a strong ev en-odd modulation, 
similar to the 1 .4 G Hz observation of [Nowakowski e t al. 
Jl982t) . At 430 MHz iBacker et alJ J 19751) found that both 
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components show an even-odd modulation, but no systematic 
drift direction was detected in the leading component. In 
our observation the 2DFS of both components of this pulsar 
contains a broad drifting feature with opposite drift sense close 
to the alias border (Fig. IA.23> . There is no evidence that the 
feature extends over the alias border, although the feature is 
not clearly separated from the alias border. 
B2021+51*: T his pulsar is clearly drifting (Fig. lA~14t . con- 
sistent with e.g. lOster et alJ il 19771) at 1720 MHz. The drifting 
is detectable over the whole P3 range. The P2 and P3 values 
that are given in tableware for the region in the 2DFS between 
the dashed lines. The centroid of the whole 2DFS gives 
P 2 = 47° + 4 and P 3 = 4.7 ± 0.1P - It is clear that the drift 
rate changes by a larg e factor during the observation, which 
wa s also observed b y lOster et all Jl977h . It was suggested 
by lOster et all i 19771) that maybe the apparent drift direction 
changes sporadically. In our observation there is no clear 
evidence that the alias mode is changing. 
B2044+15*: This observation is contaminated by interference, 
however the 2DFS of the trailing component of the pulse 
profile convincingly shows a broad drifting feature. Only the 
2DFS of the trailing component shows features and is plotted 
in Fig. |A.14I Our observation confirms the drifting found by 
Backus (1981) at 430 MHz. 



3.2.2. New drifters (classes "Dif" and "Dif*") 

B0037+56*: The 2DFS of this pulsar 5 shows a clear drift 
feature (Fig. IA.lt which appears to be extended over the alias 
border. The drift bands are visible to the eye in the pulse-stack 
and the apparent change of drift sense is also visible. There is 
also a P3 = 2fo modulation present. 

B0052+51: The trailing component of this pulsar has a broad 
drift feature in its 2DFS (Fig. IA.17l . There is a hint of drifting 
with an opposite drift sign in the first component, but this is 
not significant. The spectra also show a Pj ^ 2Pq modulation. 
B0136+57*: The drift feature is only detected in the leading 
part of the pulse profile and it appears that the feature extends 
to the horizontal axis ("Fig. IA.lt . 

B0138+59*: The drift feature in the 2DFS is broad and close 
to the horizontal axis (Fig. IA. il l. The drift feature is confirmed 
in a second observation we made. 

B0450+55*: Most of the power of the 2DFS is in the drifting 
feature (Fig. IA. 18l > and the drift bands are visible to the eye 
in the pulse-stack. The drift feature is extended to both alias 
borders. The leading component of this pulsar shows drifting 
in the opposite direction. 

B0523+11: This pulsar has a weak drift feature in the 2DFS of 
the trailing component (Fig. IA.18l . In the 2DFS of the leading 
component there is also a feature with the same P3 value, but 
in that feature there is no significant offset from the vertical 
axis measured. This means significant drifting is detected 
in the trailing component, and longitude stationary subpulse 
modulation with the same P3 value in the leadi ng compo nent. 
No drifting has been found at 430 MHz bv lBackusI Jl98ll) . but 
our observation shows that this pulsar is a drifter. 



B0525+21*: Subpulse modulation without apparent drift as 
well as some correlation between the subpulses of the two 
compone nts of the puls e profile has bee n detected for this 
pulsar bv lBackerl 1 19731) at 318 MHz and iTavlor etail ill 97-4 
at 400 MHz. We find that the two components show broad 
features to which opposite drift senses can be associated (Fig. 
IA.I8I 1. The features are possibly extended toward the P3 = 2Pq 
alias border. 

B0919+06*: The power in the 2DFS is over the whole P 3 
range is measured to be significantly offset from the vertical 
axis (Fig. IA.3t . The power in the 2DFS peaks toward the 
horizontal axis and especially this low frequency excess is 
offset from the vertical axis. This can clearly be seen in the 
bottom plot which shows a "shoulder" at the left sid e of the 
peak. No drifting has been reported for this pulsar bv lBackusI 
jl98ll) at430MHz. 

B1039-19: Both components of this pulsar show a clear, broad 
drift feature in its 2DFS with the same drift sense (Fig. IA.19l . 
B1508+55*: The subpulse modulation of this pulsar has been 
found to be unorganiz ed and without a preferred drift sense or 
a particular P3 value JTavlor & Huguenir]ll97ll at 147 MHz). 
In the 2DFS of our observation there is a broad drift feature 
present (Fig. IA.4> . which is offset from the vertical axis over 
the whole P3 range. 

B1604-00*: There are very broad features in both parts of the 
pulse profile and both components are drifting with the same 
drift sense (Fig. lA.20"l i. 

B1738-08*: The 2DFS of both halves of the pulse profile have 
a broad drift feature with the same drift sense (Fig. IA.20> . 
In the trailing component there is maybe also a second weak 
drift feature present with P3 ^ 2Pq and the same drift sense. 
The average P3 values appear to be significantly different in 
the two components, which could be because of different drift 
mode changes in the two components. The drifting can be seen 
by eye in the pulse-stack. 

B1753+52: The trailing part of the pulse profile shows a broad 
drift feature in its 2DFS (second 2DFS in Fig. IA.2U and the 
rest of the pulse profile (first 2DFS) probably has the same 
drift sense. 

B1819-22: The 2DFS of this pulsar very clearly shows a drift 
feature (Fig. IA.7L which is broadened by mode changes. A 
part of the pulse-stack is shown in Fig. ^ A full single pulse 
analysis will follow in a later paper. 

J1830-1135*: The 2DFS of this pulsar with a very long pulse 
period (6.2 seconds) shows a drift feature at the P3 = 2Pq alias 
border at +100 cpp, which is possibly double peaked (Fig. 
lA~8t . 

B1857-26: The components at both edges of the pulse profile 
are drifting with the same drift sense, which can be seen in 
Fig. IA.21I as an excess of power in the 2DFS at positive P2 
values. The drift bands are sometimes visible to the eye in the 
pulse-stack. The center part of the pulse profile does not show 
drifting in its 2DFS and is theref ore no t plotted. Th is pulsar 
is known to be a nuller jRitchingslll 97^ Biggsl fl992i) . but no 
drifting is reported in the literature. 

B1900+01*: Drifting is clearly seen over the whole P3 range 
of the 2DFS and the top part of the 2DFS is double peaked 
(Fig. IA.11> . The alias mode of this pulsar probably changes 
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during the observation. 

B1911-04*: The low frequency modulation, which is gener- 
ated in the trailing part of the pulse profile, is double peaked 
(Fig. IA. 12b . This could indicate the presence of a subpulse 
phase jump or swing or that the drift sense changes constantly 
during the observation. There seems to exist a preferred drift 
sense. 

B1917+00*: This pulsar shows a broad drifting component 
in its 2DFS, which is visible in the bottom plot as a n excess 
of pow er at positive P 2 (Fig. IA. 121 . According to IRankinI 
Jl986l) a much longer P3 50 value without a measured 
P2 was reported in a preprint by L.A. Nowakowski and T.H. 
Hankins, but to the best of our knowledge the paper was never 
published. 

B1952+29*: The drifting in both components of this pulsar is 
clearly visible to the eye in the pulse-stack and in the 2DFS 
(Fig. IA.22l . The drift sense is the same for both components. 
B1953+50*: This pulsar shows a very clear broad drifting 
feature in its 2DFS (Fig. IA.13l at low frequencies (right peak 
at P 2 = 70 cpp). 

B2053+36*: This pulsar has a broad drift feature in its 2DFS 
at low frequencies which is double peaked (Fig. IA.14I . This 
could indicate that the drift sense is changing constantly during 
the observation, which is supported by the fact that the feature 
is extended towards zero frequencies. However also a subpulse 
phase jump or swing could produce the double peaked feature. 
Subpulse modulation withou t a drift sense has been reported 
for this pulsar at 430 MHz bv lBackusI Jl98ll) . 
B2110+27*: This pulsar shows drifting over the whole P3 
range in its 2DFS (Fig. IA.15t . The lower part of the 2DFS 
is clearly double peaked, which could suggest that the alias 
mode is constantly changing during the observation. The upper 
part of the 2DFS is not convincingly double peaked. In the 
pulse-stack drifting is visible to the eye. The drift bands are 
probably distorted by nulls, causing the drift feature in the 
2DFS to be extended over the whole P3 range. Short sequences 
of drift bands with negative drifting and a P3 ^ 6P0 can be 
seen in the pulse-stack as well as some single drift bands with 
an opposite drift. A few apparent drift reversals are visible in 
the sequence, although nulling makes it difficult to identify 
them. 

B2111+46: It has been reported that this pulsar shows sub- 
pulse drift with a p ositive and negative drift sense, but without 
either dominating jTavloret alJll975l at 400 MHz). We also 
see subpulse modulation over the whole P3 range without a 
preferred drift direction in the middle and trailing components, 
but there is some systematic drift in the leading component of 
this pulsar (Fig. IA.23l . 

B2148+63*: The 2DFS of this pulsar shows broad, triple, weU 
separated features (Fig. IA. 1 5I >. The values of P2 and P3 in 
table 13 are for the feature as a whole. The centroids of the 
individual peaks give P 2 = -8.0° ± 0.2, P 3 = 2.4 ± 03P 
and P 2 = 8.3° ± 0.3, P 3 = 2.2 + 0.2P . The most likely 
interpretation of the 2DFS is that the apparent drift direction 
is constantly changing via its P3 = 2Pq alias border (see Fig. 
[2] for the expected 2DFS in this scenario). All other interpre- 
tations seems unlikely, because the feature is clearly extended 
toward the P3 = 2Pq alias border, both sides of the feature are 



separated from the vertical axis, this separation is the same on 
both sides and one side of the feature contains more power. 
The latter indicates that negative drifting dominates in this 
pulsar. 

B2154+40*: This pulsar shows a very broad drift feature 
in the 2DFS of the leading component of the pulse profile, 
but no significant drift is detected in the trailing component 
(Fig. IA.23l . The feature is probably extended toward the alias 
border. 

B2255+58: The v ery clear drift feature in the 2DFS of this 
pulsar (Fig. IA. 151 shows horizontal structure (like observed 
for instance for PSR B0809+74 and PSR B0320+39). From 
the modulation phase profile it follows that the drift bands 
make a subpulse phase step of about 140° in the middle of the 
pulse profile. The longitude resolved modulation index shows 
a minimum at the position of the subpulse phase step, as is also 
observed for PSR B0809+74 and PSR B0320+39. The phase 
step in the drift bands can be seen by eye in the pulse-stack. 
B2324+60*: This pulsar shows a broad drift feature in its 
2DFS at the alias border (Fig. IA.16l and some drift bands can 
be seen in the pulse-stack. There is also a strong P3 =s 200Po 
feature detected. 

J2346-0609*: The 2DFS of the trailing component of the 
pulse profile has a drift feature close to the alias border (Fig. 
IA.24t . The feature is not clear enough to state if the drift 
feature crosses the alias border during the observation. The 
spectra also show some low frequency modulation, especially 
in the leading component. 

B2351+61*: The 2DFS of this pulsar is double peaked at low 
frequencies (Fig. IA. 16l >. which could indicate that the drift 
direction may constantly change during the observation. Also 
the presence of a subpulse phase step or swing could produce 
this feature. The centroid is significantly offset from zero, so 
there exists a preferred drift sense during the observation. 



3.3. Longitude stationary drifters (class "Lon") 

B0402+61: The 2DFS of the trailing component of this pulsar 
shows a broad feature without a preferred drift direction (Fig. 
IA.I8I 1. The 2DFS and LRFS of the leading component is 
featureless. 

B1846-06: The 2DFS of this pulsar shows a broad feature 
with a positive value for P 2 (Fig. IA.10l i. The same drift sense 
seems to be detected at low frequencies in the 2DFS. 
B1937-26: There is no significant preferred drift sense 
detected in the 2DFS of this pulsar (Fig. IA.13l . but there seems 
to be a broad double peaked feature at the alias border. 
B1946+35: The LRFS and 2DFS shows a strong low fre- 
quency feature in both components of this pulsar (Fig. IA.22l . 
No significant offset from the vertical axis has been detected 
in the 2DFS. 

B2011+38: The broad feature in the 2DFS of this pulsar may 
have a preferred negative value for P 2 (Fig. IA.14l . The 2DFS 
and LRFS of this pulsar increases towards low frequencies and 
peaks at P 3 = 30 ± 15P - 

B2106+44: The 2DFS and LRFS of this pulsar peaks towards 
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low frequencies and seems to have a positive value for P2 (Fig. 



3.4. Unconfirmed known drifters 

B0540+23: Sporadic bu rsts of both po sitive and negative drift 
have been reported by lAshworfhl Jl982h at 400 MHz and by 
iNowakowskJll 1 99 lh at 430 MHz. No preferred drift direction 
is detected in the 2DFS of this pulsar (Fig. IA.2> . but in the 
pulse-stack short drift bands are seen with different drift senses 
confirming the previous reported drifting. Because this pulsar 
does not show a preferred drift direction in our observation, this 
pulsar is not classified as a drifter in our paper. 
B0611+22: Our observation does not show any features in the 
2DFS of this pulsar (Fig. | A.2I >. something that has also been 
reported by iBacker et alJ (1975) at 43 MHz. Drifting with 
Pt, = 50 - lOOPo has been reported bv lFerguson & Boriakoff 
( 1980) at 430 MHz, who have analyzed successive integrated 
pulse profiles. It is not clear if this kind of drifting is related 
to subpulse drifting, because in their analysis the subpulses are 
not directly measured. 

B0656+14: The modulation index of this pulsar shows a sharp 
peak at the leading edge of the pulse profile (Fig. lA.3l >. which 
is caused by a very bright subpulse. A full investigation of this 
phenomenon will be published in a upcoming paper. A pre- 
ferred negative drift sense has been reported by Backus: ( 1981) 
at 430 MHz. We do not see a preferred drift sense, but there is 
low frequency modulation. 

B0820+02: Positive drifting has been reported by iBackusj 
Jl98ll) at 430 MHz. Our observation probably lacks the S/N 
to confirm this IA. 31 

B0950+08: Drifti ng has been re ported for this pulsar (e.g. 
iBackeJI 19731 and IWolszczanll 1 980) with a variable P 3 - 6.5. 
This drifting is not visible in the 2DFS of our observation (Fig. 
IA.4> . but subpulse modulation is seen over the whole P 3 range 
without a preferred drift sense. The interpulse has no mea- 
su red mo dulation and is not plotted. The observing frequency 
of Backer) Jl973l) was 430 MHz, so it could be that the drifting 
of this pulsar is only vi sible at low observ ing frequencies. The 
observing frequency of lWolszczanl Jl98d) is not mentioned in 
their paper. 

B1112+50: This pulsar is known to show nulling and pulse 
profile mode switching an d in one of these modes drifting 
subpulses are reported (e.g. Wri ght et alJI 19861 at 1412 MHz). 
There is no clear drift feature detected in the 2DFS of our ob- 
servation, but subpulse modulation is seen over the whole P 3 
range (Fig. IA.19t . 

B1612+0 7: N egative subpulse drift has been reported by 
iBackusI lll98lh at 430 MHz for this pulsar 5 . The 2DFS of our 
observation is featureless (Fig. IA.4l . which could be because a 
too low S/N. 

B1918+19: This pulsar is show n to be a drifter with at leas t 
four drift modes at 430 MHz (Hankins & Wolszczanl fl987). 
There are no features in the 2DFS of our observation (Fig. 
IA. 12i . which could be because a too low S/N. 
B2315+21: Drifting with a negative drift sense has been re- 



ported for this pulsar at 430 MHz bv lBackusNl98ll) .Our spec- 
tra (Fig. IA.16l l do not show any sign of drifting, what could be 
because a too low S/N. 

3.5. Pulsars with low frequency modulation 

A modulation index could be derived from the observations of 
the following pulsars and their spectra show an excess of power 
toward the horizontal axis (i.e. a "red" feature). This means 
that there is some correlation between successive pulses, but 
no quasiperiodicity. 
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B1804-08: The 2DFS of this pulsar possably shows a broad 
drift feature which is generated primarily by the trailing com- 
ponent of the pulse profile ("Fig. IA.6t . but after randomizing the 
order of the pulses this feature turned out to be not significant. 
B1924+16: The 2DFS of this pulsar shows a hint of a broad 
drifting feature (Fig. IA. 1 31 . but scrambling the pulse stack 
showed that this drifting is not significantly detected. 



3.6. Pulsars with a flat spectrum 

A modulation index could be derived from the observations of 
the following pulsars and their spectra show no clear features 
over the vertical range. This means that the subpulse modula- 
tion appears disordered or that the S/N ratio is too low. 



B0144+59: Two vertical bands are detected in the 2DFS of 
the trailing component (best visible in the bottom plot of Fig. 
IA.17I . This bands may also be present (although weakly) in 
the middle component. A Pi — ±2.3° (160 cpp) subpulse 
separation can be associated with this feature, but no particular 
P3 value. This indicates that there is a quasiperiodic intensity 
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modulation in the pulses with a period of about 1.3 ms, 
but there is no correlation in the positions of the subpulses 
from pulse to pulse. The same features are seen in another 
observation we made of this pulsar. We see the same kind 
of phenomenon (a bit more clear) for PSR B2217+47. The 
leading component does not show any features and is therefore 
not plotted. The spectra shown are calculated using transforms 
of only 32 pulses in order to reduce the resolution. This makes 
it more easy to see the features by eye. 

B0531+21: Both the 2DFS of the main and interpulse of the 
Crab pulsar does not show any sign of drifting (Fig. lA.1%1 . A 
very large modulation index measured is measured (m _ 5), 
which is caused by its giant pulses JStaelin & ReifensteirJ 
1968). 

B0626+24: The 2DFS of this pulsar shows subpulse modula- 
tion over the whole P3 range (Fig. IA.2t . Subpulse modulation 
without a drift sense has been reported by iBackusI dl98lh at 
430 MHz. 

J1022+1001: This millise cond pulsar is know n to show 
subpulse modulation (Edwards & Stappers 2003a). The power 
in the 2DFS of the trailing component pea ks toward P3 = 2Pp 
(Fig. | A.4> . consistent with the analysis of Edwards & Stappers 
1 2003a) at 1380 MHz (we have used the same data). The 2DFS 
shown in the figure is that of the trailing component. 
B1541+09: This pulsar (Fig. IA.4> is observed to have a low 
frequency excess and exhibits mo de changes and orga nized, 
but short, drifts in both directions (Nowa kowskill 19961 at 430 
MHz). 

J1713+0747: This pulsar shows subpulse modulation 
(Edwards & Stappers1 l2003al at 1190 and 1700 MHz), but the 
quality of our observation is too poor to confirm this (Fig. 

ii~5i . 

B1736-29: The interpulse of this pulsars is not plotted, be- 
cause no features are detected in its spectra and no modulation 
index has been measured (Fig. lA.5> . 

B1737+13: This pulsar shows a clear P3=ll-14Po longitude 
stationary subpulse modulation, but no drifting has been 
detected iRankin et al.ll988l at 1412 MHz). No features appear 
in the spectra of our observation (Fig. |A.5I . No drifting has 
been detected for this pulsar by BackusJ Jl98ll) at 430 MHz. 
B1749-28: This pulsar shows flat featureless spectra (Fig. 
IA.6i . A flat featureless fluctuation spectrum has also been 
observed for this pulsar at a lower observing frequency by 
iTavlor & Hugueninll97ll 

B1818-04: The power in the 2DFS is possibly double peaked 
(Fig. IA.7l . This could be because of the presence of a subpulse 
phase step or swing or because of drift reversals. It has been 
reported that the subpulse modulatio n is not we ll organized 
jTavlor & Hugueninl 197 H and ITavlor et all i 19751) both at 400 
MHz). 

B1839+09: Subpulse modulation without any drift sense has 
been detected by Backus (1981) at 430 MHz. No features 
appear in the spectra of our observation (Fig. lA.9> . 
B1842-04: When we first analyzed the spectra of this pulsar, 
an extremely bright and surprisingly sharp longitude stationary 
P3 = 3.00Po subpulse modulation feature appeared. Folding 
the data with a three times longer pulse per iod revealed that th e 
pulse period of this pulsar, as reported by Clift on et alJ (|l,992), 
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Fig. 4. The fraction of pulsars we observe to show the drift- 
ing phenomenon (solid line) and the number of pulsars (dashed 
line) versus the measured S/N ratio of the observation. The 
root-mean-square (RMS) is calculated as an estimate for the 
error (if the bin contains more than one observation). The dot- 
ted line is a fit for the S/N dependence of the chance to detect 
drifting subpulses. 



is wrong by a factor three. It turns out that the correct pulse 
perio d of this pulsar has appeared in the literature jHobbs et alJ 
2004) without a comment about this discrepancy. Private 
communication with G. Hobbs revealed that this deviation 
is first discovered, although apparently not reported in the 
literature, by the Parkes multibeam survey. Using the correct 
pulse period the spectra are featureless (Fig. IA.10l . 
B1842+14: Subp ulse modulatio n without a drift sense has 
been detected bv lBackusI Jl98lh at 430 MHz. The spectra of 
our observation is featureless (Fig. IA.10l . 
J1850+0026: The shown 2DFS of this pulsar (Fig. lA~TO is of 
the trailing peak. 

B2022+50: The interpulse of this pulsar (which also does not 
show any features in its spectum) is not plotted (Fig. IA.14l . 
J2145-0750: Weak quasi-periodicities around 0.22 and 
0.45 cpp are visible in the LRFS and 2DFS of the leading 
component of this milli second pul s ar (Fig . IA. 15l >. as has been 
reported by lEdwards & Stappers! d2003al) at 860 and 1380 
MH z. We have used the sam e 21 cm data as has been used by 
Edwards & Stappers (2003a). 

B2217+47: There is no preferred drift sense detected in the 
feature in the 2DFS of this pulsar (Fig. IA.15II. which would 
confirm the observation of Taylor & Hugueninl 1 1 97 ll at 147 
MHz. However the 2DFS shows two vertical bands smeared 
over the whole Pj range. This modulation is primarily gener- 
ated in the right part of the pulse profile. The interpretation is, 
like for PSR B0 144+59, that there is a quasiperiodic intensity 
modulation in the pulses with a period of about 2.5 ms, but 
there is no correlation in the positions of the subpulses from 
pulse to pulse. 
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Fig. 5. The P-P diagram of the analyzed pulsars (including the low S/N observations), where P is the pulse period and P its 
time derivative. The non drifting pulsars are the dots, the diffuse (Dif and Dif) drifters are the open stars, the coherent drifters 
are the filled stars and the pulsars showing longitude stationary subpulse modulatio n are the open circles. Line s of equal surface 
magnetic field strength and characteristic ages are plotted, as well as a death line lichen & Rudermanl ll993h . The millisecond 
pulsars are not plotted to make the plot more readable. 



3.7. Pulsars without a measured modulation index 

No modulation index could be measured from our observations 
of these pulsars, so no subpulse modulation could be detected. 
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B1821+05: Subpu lse modulati on without a drift sense has 
been reported by B ackus! dl98ll) at 430 MHz. No features are 
seen in the 2DFS (Fig. IA.8t and no modulation index could be 
measured for this pulsar. This is probably because of the low 
S/N of our observation. 

B1915+13: No feature s are seen in the spectra of this pulsar 
by backer et al.l (119751) at 430 MHz. In our observation there 
are also no features (Fig. IA.12l . which could be because of the 
low S/N of our observation. 

B1937+21: The spectra of this pulsar are featureless and there 
is no modulation index measured. Only the 2DFS of the main 
pulse is plotted in Fig. IA.13l 



4. Statistics 

4.1. The numbers 

The selection of our sample of pulsars is based only on the pre- 
dicted S/N in a reasonable observing time. While this sample is 
luminosity biased, it is not biased on pulsar type or any partic- 
ular pulsar characteristics. This allows us, first of all, to address 
the very basic question: what fraction of the pulsars show the 
drifting phenomenon? 

Of the 187 analyzed pulsars 68 pulsars show the drifting 
phenomenon, indicating that at least one in three pulsars drift. 
This is however a lower limit for a number of reasons. First 
of all, not all the observations have the expected S/N. This 
could be because of interference during the observation, inter- 
stellar scintillation, digitization effects, or because the flux or 
pulse width for some pulsars was wrong in the database used. 
The latter was confirmed by measuring the pulse with directly 
from our own observation. Also there are 6 pulsars which show 
longitude stationary subpulse modulation. Longitude station- 
ary subpulse modulation could indicate that there is drifting, 
but without a preferred drift sense and therefore it could be re- 
lated to the same phenomenon. 

Because many pulsars in our sample were found to be drift- 
ing, the sensitivity of our method to detect drifting could be 
checked. The S/N of the observations is determined by com- 
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Fig. 6. The left panel shows the histogram of the characteristic ages of the analyzed pulsars with a S/N > 100. The solid line 
is the age distribution of the non-drifting pulsars, the dashed line shows all the drifters and the dotted line shows the coherent 
drifters. The right panel shows the "S/N versus age bias" corrected histogram. 



paring the root-mean-square (RMS) of the off-pulse region of 
the pulse profile with the power in the pulse. The width of the 
pulse was automatically determined by trying different pulse- 
widths and maximizing the resulting S/N. In Fig.0]the fraction 
of pulsars that show the drifting phenomenon is plotted versus 
the S/N ratio of the observation. One can see that the probabil- 
ity of detecting drifting is higher for observations with a higher 
S/N. The method used is also working for observations with a 
low S/N, but not in all cases. This could be because, in order 
to detect drifting in observations with a low S/N, the fraction 
of pulse energy that is in the drifting subpulses should be high. 
Also the fraction of time the pulsar shows drifting subpulses 
during the observation should be high and the drifting must be 
reasonably coherent. 

To make the statistics more independent of the S/N ratio of 
the observations, we will do the statistics for the observations 
with a S/N > 100. For 106 pulsars in our sample this criterion 
is met. By excluding observations with a low S/N a bias toward 
long period pulsars (which are observed longer to get enough 
pulses) and well studied pulsars (for which long archival data 
was available) may be introduced. Therefore all the statistics 
are checked including the low S/N observations. 

Besides the S/N of the observations, the amount of scatter- 
broadening could also influence the probability to detect drift- 
ing subpulses. If the amount of scatter broadening is wider than 
the subpulse separation P2, then the sensitivity to detect drift- 
ing subpulses will be severely diminished. There are a number 
of pulsars in our sample that seem to show scatter broadening. 
It is however difficult to distinguish between a pulse profile 
that is scatter broadened and a profile which has an intrinsic 
exponential tail like shape. Most pulsars that probably show 
scatter broadening have a S/N below the threshold value of 
100, so they are therefore excluded from the statistical anal- 
ysis (PSRs B1758-23, B1817-13, B1822-14, J1828-1101, 
B 1832-06 and B 1849+00). There is only one pulsar that seems 



to show scatter broadening and has a S/N above the threshold 
value (PSR B1815-14). Because the low number of pulsars 
that are scatter broadened, it seems very likely that their influ- 
ence on the statistics can be neglected. Scatter broadening will 
be more of an issue for our subsequent paper, which will focus 
on a lower frequency study of the subpulse modulation proper- 
ties of radio pulsars. 

Of the pulsars with high enough S/N observations, 57 are 
detected to be drifters (54%) and 5 pulsars show longitude sta- 
tionary subpulse modulation (5%). From Fig. 0] it is clear that 
the real drift percentage could even be hig her than 54%. Thi s 
number is consistent with As hworthl l l 19821) and lBackusMl98ll) . 
who found about the same number based on a smaller sample 
of pulsars. There are many reasons why drifting is not expected 
to be detected for all pulsars. For instance for some pulsars 
the line of sight cuts the magnetic pole centrally and therefore 
longitude stationary subpulse modulation is expected. Also, re- 
fractive distortion in the pulsar magnetosphere (e .g. IPetroval 
2000| IWeltevrede et~ai]l2003l iFussell & Luoll2004 or nulling 
will disrupt the drift bands, making it difficult or even impos- 
sibl e to detect drifting. A = 50 - 100Po has been reported 
bv lFerguson& Boriakofj|ll980l) for PSR B0611+22, indicat- 
ing that the P3 value for some pulsars could be very large. In 
that case longer observations are needed to detect this drifting 
and distinguishing it from interstellar scintillation could be- 
come a problem. Some pulsars are known to show organized 
drifting subpulses in bursts. In that case some of our observa- 
tions could be too short to contain enough drift bands to detect 
the drifting. 

With a lower limit of one in two it is clear that drifting is 
at the very least a common phenomenon for radio pulsars. This 
implies that the physical conditions required for the emission 
mechanism of radio pulsars to work cannot be very different 
than the physical conditions required for the drifting mecha- 
nism. Therefore it could well be that the drifting phenomenon 



P. Weltevrede et al.: The subpulse modulation properties of pulsars at 21 cm 



19 



is an intrinsic property of the emission mechanism, although 
for some pulsars it is difficult or even impossible to detect. 

4.2. The drifting phenomenon and the P-P diagram 

The unbiased sample of pulsars not only allows us to determine 
what fraction of the pulsars show the drifting phenomenon, 
but also to correlate the drifting phenomenon with other pulsar 
parameters. Two directly measurable and therefore important 
physical parameters of the pulsar are the pulse period and its 
time derivative (spin-down parameter). From the position of a 
pulsar in the P-P diagram and assuming magnetic dipole brak- 
ing, an estimate of the age and the magnetic field strength can 
be obtained. Therefore it is useful to try to correlate the drift- 
ing phenomenon with the position of the pulsars in the P-P 
diagram (Fig. [5}. All the analyzed pulsars with a measured P 
are in this diagram 8 and the coherent drifters, diffuse drifters 
and pulsars showing longitude stationary subpulse modulation 
are plotted with different symbols to identify any trends be- 
tween position and classification of the pulsars. To make the 
plot more readable the millisecond pulsars are not plotted. 

The P-P diagram reveals that the pulsars that show the 
drifting phenomenon are more likely to be found closer to the 
death line and this is more pronounced for the coherent drifters. 
This suggests that the population of pulsars that show the drift- 
ing phenomenon is on average older (the pulsar age is defined 
as t — jP/P) than the population of pulsars that do not show 
drifting. This confirms the result of lAshworthl i 1982b . who also 
found that drifters are on average older. Moreover it seems that 
drifting is more coherent for older pulsars. This trend is more 
visible in the pulsar age histograms (left panel of Fig.|6j, where 
the nondrifters, the coherent drifters and all the drifters (both 
coherent and diffuse) are plotted separately. This correlation 
seems to suggest an evolutionary trend that the subpulse mod- 
ulation is disordered for the youngest pulsars and gets more and 
more organized into drifting subpulses as the pulsar ages. 

The significance of this trend can be determined with the 
Kolmogorov-Smirnov test (KS-test), which tells us how likely 
it is that two distributions are statistically different. It follows 
that the age distribution of the drifters is only 0.03% likely to be 
the same as the age distribution of pulsars not showing the drift- 
ing phenomenon. Thus the drifters and nondrifters have signif- 
icantly different age distributions. The KS-test is also used to 
find out if the coherent drifters have a separate age distribu- 
tion. It follows that the coherent drifters are only 0.4% and 8% 
likely to have the same age distribution as the nondrifting and 
the diffuse drifting pulsars respectively. Therefore the pulsars 
which drift coherently are likely to have a separate age distri- 
bution. Although likely, the difference in the age distribution of 
the coherent drifters is not detected to be significantly different 
from the drifters. Nevertheless it is intriguing to think that drift- 
ing becomes more and more coherent for pulsars with a higher 
age. A larger sample of pulsars is needed to check whether this 
is significant. The same trend is found when the low S/N ob- 
servations are included in our sample. In that case the drifters 
and coherent drifters are respectively 0.03% and 0.1% likely 
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Except PSR B1744-24A, which has a negative P. 



Fig. 7. The surface magnetic field strength histogram of the 
pulsars which do not show the drifting phenomenon (solid 
line), those which do show the drifting phenomenon (dashed 
line) and those which drift coherently (dotted line). Only ob- 
servations with a S/N > 100 are included. 



to have the same age distribution as the nondrifters, confirming 
the assertion that the drifters and non-drifters have different age 
distributions. 

The S/N of the observations used in the left panel of Fig. 
[6] may be different in each of the age bins, thereby introduc- 
ing a "S/N versus age" bias. To correct for this effect, for each 
distribution and for each age bin the median of the S/N of the 
observations was calculated. This median S/N was used to es- 
timate what the chance was of detecting drifting (the dotted 
line in Fig. 0}. The distributions of the drifters and coherent 
drifters were divided by this chance (because a low S/N implies 
that there was only a low chance of detecting drifting) and the 
distribution of the non drifters was divided by one minus this 
chance (because a low S/N implies a high chance of not detect- 
ing drifting). This gives the "S/N versus age bias" corrected 
age distributions (right panel of Fig. [6jl. As one can see this 
correction does not lead to a qualitatively different result. 

A possible explanation for the age dependence of the drift- 
ing phenomenon is that the drift bands are more distorted for 
younger pulsar s. One m echanism to dis t ort the drift b ands 
is nu lling (e.g.lTavlor&Hugueninlll97ll Ivan Leeuwen et alJ 
l2002t Ijanssen & van Leeuwenll2004l) and the fraction of time 
that pulsars spend in their nulling state (the nulling fraction) is 
known to be correlated with the pulsar age. However it has been 
found bv lRitching sl i 197 (h that the nulling fraction is on aver- 
age higher for olde r pulsars, w hich is confirmed in later studies 
(e.g. Li & Wang 1995; Biggs 1992). Although the correlation 
with other pulsar parameters seems to be stronger, these studies 
prove that nulling cannot explain this correlation. 

Another possible scenario is that the alignment of the mag- 
netic dipole axis with the rotation axis has something to do with 
the observed trend. Observations seem to show that the angle 
a between the magnetic axis and the rotation axis is on aver- 
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age smaller for older pulsars (e.g. lTauris & ManchesteJ 1*9 98). 
indicating that the magnetic axis and the rotation axis becomes 
more aligned or anti-aligned for older pulsars. This angle is 
likely to be an important physical parameter in the mechanism 
that driv es the drifting phenomen on (for instance the classi- 
cal iRuderman & Sutherlandin"975l model can only be applied 
for an anti-parallel magnetic axis). In this scenario as the pul- 
sar gets older, the rotation axis and the magnetic axis grows 
more aligned, which makes the drifting mechanism more ef- 
fective or regular. This trend is also consistent with the fact 
that pulsar s with a regular drift pattern te nd to have small val- 
ues for a dWrightl [20031 lRankiniri993bl) . However, we have 
found that the interpulse pulsar PSR B 1702- 19 is a coherent 
drifter, suggesting that coherent drifters ca n have a large a (an 
a value of 85° and 90° hav e been found by Rankin 1993bjand 
iLvne & ManchesteJ 1 19881 respectively). Also t he pulse mor- 
phology seems to e volve when the pulsar ages jRankinll 19831 
Besk inet al]ll984 . such that core single stars are on aver- 
age younger than pulsars with more complex profiles. This 
could make drifting subpulses more likely to be detected in 
older pulsars. In the non-radial pulsations model this trend can 
also be explained, because the appearance of narrow drifting 
subpulses is favored in pulsars with an aligned magnetic axis 
jClemens & Rosenl2004 . 

An estimate for the component of the surface magnetic field 
of pulsars perpendicular to the rotation axis can be directly 
derived from the position of the pulsar in the P-P diagram 
(B s = lO^VTo 15 PP Gauss). The histograms of the magnetic 
field strengths of the three different groups of pulsars (Fig.0 
do not show a clear trend, which is confirmed by the KS-test. 
It follows that the magnetic field strength distribution of the 
nondrifters has a chance of ~20% and ~50% to be statistically 
the same as the distributions of all the drifters (both the coher- 
ent and diffuse drifters) and the coherent drifters respectively. 
This means that the magnetic field strength distributions are not 
significantly different. If the low S/N pulsars are included, the 
magnetic field distributions are more likely to be the same than 
to be different. 

It seems that the drifting phenomenon is only weakly cor- 
related with, or even independent of magnetic field strength. 
This is consistent with the large fraction of pulsars that are 
found to show the drifting phenomenon, because the drifting 
phenomenon is too common to require very special physical 
conditions. 



4.3. The drifting phenomenon and the modulation 
index 

The drifting phenomenon is a form of subpulse modulation, so 
the longitude-resolved modulation index m, (Eq.[3) is an obvi- 
ous parameter to try to correlate with the drifting phenomenon. 
Because the longitude-resolved modulation index can vary a lot 
with pulse longitude, as can be seen in the figures in appendix 
[A] it is a somewhat arbitrary what one should call the modula- 
tion index. The longitude-resolved modulation index of many 
pulsars do show a minimum in the middle of the pulse pro- 
file where the total intensity is relatively high. This means that 
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Fig. 8. The modulation index distribution of the pulsars that 
do not show the drifting phenomenon (solid line), that do show 
the drifting phenomenon (dashed line) and of the pulsars that 
drift coherently (dotted line). Observations with a S/N < 100 
and PSR B0531+21 (with a measured m — 5) are not included 
this plot. 



if the S/N of an observation increases, the average modulation 
index will also increase. This is because then the modulation in- 
dex can be measured at pulse longitudes farther away from the 
peak intensity of the pulse profile where the longitude-resolved 
modulation index tends to be higher. We have therefore chosen 
the modulation index m to be the longitude-resolved modula- 
tion index m, at the pulse longitude bin i where m, has its mini- 
mum value. This definition should make the modulation index 
a more S/N independent number than for instance the average 
of the longitud e-resolved m odulation index. The same defini- 
tion is used by lTenetfeGill ll2003) to measure the modulation 
index. 

Modulation index histograms are shown in Fig. [8] Readily 
apparent is the trend that pulsars that show the drifting phe- 
nomenon more coherently have on average a lower modulation 
index. There seems to be no significant difference in the mod- 
ulation index of the pulsars that do and do not show the drift- 
ing phenomenon. The significance of this trend is also checked 
with the KS-test and the modulation index distribution of the 
drifters is ~5Q% as likely to be the same as the distribution 
of the nondrifting pulsars. The coherent drifting distribution is 
only 6% and 5% likely to be the same as the nondrifters and 
the diffuse drifting distributions respectively. These numbers 
are too high to state that the modulation index distribution of 
the coherent drifters are significantly different. 

While the trend seen in Fig. [8] is not shown to be statisti- 
cally significant, it is intriguing and a larger sample of pulsars 
is needed to determine its true significance. If the correlation is 
proven to be significant it would indicate that pulsars that show 
coherently drifting subpulses have on average a lower modu- 
lation index. Although this trend may appear counterintuitive 
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because drifting subpulses imply subpulse modulation, it can 
be explained. 

From the summation in Eq.0](and demonstrated in Fig.[3J 
it follows that the modulation index is independent of whether 
the subpulses appear randomly or organized and from Eq. [5] 
it follows that the modulation index is even independent of 
the drift band separation P3 in the case of a coherent drifter. 
However, if the number of subpulses per pulse is large, the sub- 
pulses could possibly overlap causing the intensity to change 
less from pulse to pulse resulting in a lower modulation index 
ilenet&Gi]l E003). It must be noted that since sparks cannot 
physically overlap on the polar cap, this only works if there is 
significant broadening in the mapping from polar cap to the ra- 
diation beam pattern. Another parameter influencing the mod- 
ulation index is the width of the subpulse intensity distribution. 
If the subpulses have a narrow subpulse intensity distribution, 
e.g. the subpulses have more equal intensities, the modulation 
index will also be lower. A clear example of this effect is the 
huge measured modulatio n index of P SR B053 1+21 ( m = 5), 
caused by its giant pulses llStaelin & Reifen stein 196a). 

To explain the trend that pulsars that show the drifting 
phenomenon coherently have on average a lower modulation 
index, these pulsars must either have on average more sub- 
pulses per pulse or the subpulse intensity distribution must be 
more narrow. In the sparking gap model it seems reasonable 
that the subpulses of the coherent drifters have more equal 
intensities. Coherent drifting could indicate that the electro- 
dynamical conditions in the sparking gap are stable, which 
could be the reason why the subpulses have on average more 
equal intensities. Also the presence of subpulse phase steps re- 
sults injyninirrmm injh^ in- 
dex JEdwards et aljEool lEdwards & Stannersll2003bl). This 
effect can be seen in the longitude resolved modulation index of 
PSR B0320+39, PSR B0809+74, PSR B1919+21 and the new 
drifter PSR B2255+58. Also PSR B0818-13 shows a mini- 
mum in its longitude resolved modulation ind ex at the posi- 
tion of its subpulse phase swing. It is argued by Ed wards et alJ 
(2003) that the local reduction of the modulation index accom- 
panied by a rapid swing in the modulation phase profile are the 
result of interference between two superposed drifting subpulse 
signals that are out of phase. It is not unlikely that interference 
can only occur if the drifting is coherent, which could explained 
the trend. 

It should also be noted that the modulation index of a purely 
sinusoidal subpulse signal results in a modulation index of 
1 / V2. Subpulse patterns with different waveforms or drift band 
shapes will generally have larger modulation indices. Many 
pulsars have a modulation index which is significantly lower 
than this value. This implies that the pulsar emission has both 
a subpulse signal and a non-varying component, which could 
indicate the presence of superposed out of phase subpulse sig- 
nals. 

Another explanation for this trend would be that refraction 
is perhaps more dominant for pulsars that do not show coher- 
ently drifting subpulses. The pulse morphology could well be 
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Fig. 9. The average modulation index versus pulsar age 
histogram for the pulsars showing the drifting phenomenon 
(dashed line) and all the pulsars with a measured modulation 
index (solid line). The RMS is calculated as an estimate for the 
error (if the bin contains more than one pulsar). Pulsars with 
a S/N < 100 and PSR B0531+21 (r = 1240, m = 5) are not 
included in this plot. 



the organized drifting subpulses are more refractively distorted 
than for others. For those pulsars the coherent drifting is dis- 
torted in this scenario, causing the subpulses to appear more 
disordered in the pulse window. Moreover it is expected that 
the intensities o f the individua l subpulses varies mor e because 
of lensing (e.g. IPetrovalboOOt iFussell & Luoll2004) and pos - 
sible focusing of the radio emission ( Weltevred e et alJ l2003t) . 
causing the modulation index to be higher in those pulsars. 

If a correlation between the drifting phenomenon and both 
the pulsar age and the modulation index exist, there could also 
be a correlation between pulsar age and modulation index. The 
modulation index versus pulsar age histogram is plotted in Fig. 
[9] and it is clear that no significant correlation is found, indi- 
cating that the modulation index is the same for pulsars with 
different ages. This seems to suggest that a high pulsar age and 
a low modulation index are two independent factors affecting 
the likelihood that a pulsar will exhibit coherently drifting sub- 
pulses. 

The m odulatio n inde x of core type emission is observed 
by IWeisberg et alJ Il986l to be in general lower than that of 
conal type of emis sion. This is also a consequence of the 
Gil & Sendvk 2000 model. In the sparking gap model, the drift- 
ing phenomenon is associated with conal emission and there- 
fore expected to be seen in pulsars with an on average higher 
modulation index. The modulation index distributions of the 
drifters and nondrifters are likely to be the same, so drift- 
ing phenomenon appears not to be an e xclusively conal phe- 
nomenon (as suggested by Rankin 1986). If well organized co- 
herent drifting is an exclusively conal phenomenon, it is ex- 
pected that coherent drifters have an on average a higher mod- 
ulation index, exactly opposite to the observed trend. 
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Fig. 10. The modulation index for all analyzed pulsars which have a measured modulation index (except PSR B0531+21) versus 
the four complexity parameters as described in the text. 



Also, if drifting is an exclusively conal phenomenon, no 
drifting is expected for pulsars classified as "core single stars". 
Although this may be true for many cases there are some excep- 
tions. The interpulse of PSR B1702-19 is classified as a core 
single star (Rankin 1990 and references therein) and shows a 
clear and narrow P3 feature. The diffuse drifter B2255+58 is 
another good example of a drifter that is classified as a core 
single star. Most of the core single stars that show drifting are 
diffuse Dif drifters (PSR B0136+57, B0823+26, B1642-03, 
B 1900+01, B1911-04, B1953+50 and B2053+36). The co- 
herent drifter PSR B 1844-04 has been classified as a core sin- 
gle or a triple profile. This means that it is questionable if the 
lack of of ordered subpulse modulation is a use ful criterion to 
identify core emission (as suggested by Rankin 1986). It must 
be noted that the classification can be frequency dependent, so 
core single pulsars at low frequencies could show conal emis- 
sion at higher frequencies. The many core single stars that ap- 
pear to be drifting stresses the importance of being unbiased on 
pulsar type when studying the drifting phenomenon. 



4.4. Complexity parameter 

In the framework of the sparking gap model (e.g. 
Rud erman & Sutherland! fl97l iGil & Sendvkll2000t iGil et alJ 
2003) the subpulses are generated (indirectly) by discharges 
in the polar gap (i.e. sparks). Each individual spark should 
emit nearly steady, unmodulated radiation, so the modulation 
of the pulsar emission is due to the changing positions of 
the subpulses in the pulse window and the number of visible 
sparks in different pulses. The more sparks there are visible in 
the pulse window, the less the intensity will change from pulse 
to pulse because the subpulses could overlap. The number of 
sparks that fits on the polar cap is quantified by the complexity 
parameter & Sendvk 2000) and therefore one expects an 
anti-correlation between the modulation index m (which is a 
measure for how much the inte nsity varies from pulse to pulse) 
and this complexity parameter djenet & Gill2003l) . As noted in 
the previous subsection, this only works if there is significant 
broadening in the mapping from polar cap to the radiation 
beam pattern. 

The complexity parameter is a function of the pulse pe- 
riod and its derivative and its precise form depends on the 
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Fig. 11. The probability functions P(p) of the correlation coef- 
ficient p for the four complexity parameters, the pulse period, 
the spin-down parameter and the age of the pulsar. The posi- 
tion of the maximum as well as the 1-cr widths are indicated as 
well. 



model one assumes for the pulsar emission. By correlating the 
modulation index of a sample of pulsars with various com- 
plexity parameters as predicted by different emission models 
one could try to distinguish which model best fits the data 
Jjenet & Gill2003l) . We have correlated the modulation indices 
in our sample of pulsars with the comp lexity parameter of four 
different emission mo dels as derived bv l.Tenet & Gill J2Q03ll and 
lGil&Sendvklll2000l) : 



fli = 5(P/10- I5 ) 2/7 (P/ls)- 9/14 , a 2 = VPP- 3 
a 3 = Vpp, fl4 = 



(7) 



These are the complexity p arameters for the sparking gap 
model (Gil&Sendvk 2000), continuous current outflow in- 



2001 



2001) 
2003 



stabilities ( Arons & Scharlemann 1979; Hibschman & Arons 
surface magnetohydrodynamic wave insta bilities Ilou 
and outer magnetospheric instabilities Jjenet & Gi 
) respectively. 



Physically, a\ is proportional to the polar cap radius r p di- 
vided by the gap height h as predicted in the sparking gap 
model and is therefore a measure for the number of sparks 
across the polar cap. The parameter a 2 is proportional to the 
acceleration parameter, which is the surface magnetic field 
strength B s divided by P 2 . This acceleration parameter is pro- 
portional to total current outflow from the polar cap, and 
roughly to the circulation time of the sparks expressed in pulse 
periods. The acceleration parameter is also proportional to the 
square root of the spin down energy loss rate. Finally 03 and a\ 
are proportional to respectively the magnetic field strength at 
the surface and at the light cylinder. 

According to Ijenet & Gill |2003) the anti-correlation be- 
tween the modulation index and the complexity parameter will 
be masked by viewing angle issues in conal emission, so one 
might get better results by including only pulsars that are 
known to emit core emission. Because, as discussed in Sect. 
14.31 it is not clear how drifting subpulses relate to the morpho- 
logical classification of the pulsar, all pulsars with a measured 
modulation index are included in our sample. The modulation 
index is chosen to be th e minimum in the longitude-resolved 
modulation index (like in iTenet & Gi]l2003l which should give 
the best estimate for the modulation index of the core emission 
if present in the pulse profile. 

The modulation index versus the four complexity param- 
eters plots are shown in Fig. ^| To find out if there exist an 
(anti-)correlation without fitting a specific function to the data 
a rank-order correlation is used. This means that the rank of the 
values among all the other values is used rather than the values 
itself. This implies that the correlation coefficient is identical 
for the set of points (xu yd and (F(xt), G(y,)), as long as the 
functions F and G are monotonic functions. This means that 
for instance, because a 2 is proportional to the square root of 
the spin down energy loss rate, the correlation coefficient of a 2 
and the modulation index will be same as the correlation coef- 
ficient of the the spin down energy loss rate and the modulation 
index. 

Following Lfenet & Gill d2003i) we have used the Spearman 
rank-order ed correlation co efficient p and its significance pa- 
rameter A ( Pre ss et al.ll9 92). A problem arises when one wants 
to include the uncertainties of the data points, because then the 
rank of the values is not uniquely defined anymore. As one can 
see in Fig. [TO] the errorbars are overlapping each other, so they 
should be included in the analysis. The significance parameter 
A does not include the uncertainties on the data points and is 
therefore not directly usable to estimate the significance of p. 

To include the uncertainties in the analysis we have used 
a Monte Carlo approach. The data points are replaced with 
Gaussian distributions with a width corresponding to the 1-cr 
uncertainties of the measurements. In each integration step a 
point is randomly picked from these distributions. Instead of 
calculating p and A directly from the data points, we calculate 
them for the randomly choosen points. So for each integration 
step we randomly select a set of points for which we get a p 
and A. The probability distribution P(p) is calculated by aver- 
aging the Gaussian distributions centered around the calculated 
values of p with a 1-cr width A. The calculated probablility dis- 
tributions are plotted in Fig.^J The position of the peak of the 
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Table 1. The correlation coefficients and their significance as 
derived from Fig. II II 



probability distribution corresponds to the most likely value of 
the correlation coefficient and the l-cr width of the peak is a 
measure for the significance of the correlation coefficient. 

The results of this analy sis are tabu l ated in tabled Based 
on a sample of 12 pulsars. Ijenet & Gill ( 120031) concluded that 
the sparking gap model (ai) showed the highest anti-correlation 
and that the surface magnetohydrodynamic wave instabilities 
(«3) is unlikely. Also in this enlarged sample, 03 shows the 
least evidence for an anti-correlation (it is even more likely 
that the modulation index is positively correlated with aj). 
The strongest anti-correlations are found for and a\, which 
corresponds respectively to continuous current outflow insta- 
bilities and outer magnetospheric instabilities. However, none 
of the correlations are significantly inconsistent with an anti- 
correlation, and therefore none of the models can be ruled out 
based on these observations. 

One can also see that the modulation index is uncorrected 
with the age of the pulsar, which is consistent with Fig. [9] There 
is a hint that the modulation index is weakly correlated with the 
pulse period and the surface magnetic field strength B s . 

4.5. Properties of drift behavior 

The value of Pt, is observed to be i ndependent of the observ- 
ing frequency Jlzyekoya $\ al]|l 993|). but the value of Po could 
vary a little (e.g. Edwards & S tappers 2003a). Moreover obser- 
vations sho w that measuring a value f or P2 can be far from 
trivial (e.g. E dwards & Stappersl2003bl) and it is only a mean- 
ingful parameter if the drift bands are linear. This means that 
correlating P3 with other pulsar parameters is the most direct 
way to find out if the drift rate depends on any physical param- 
eters of the pulsar. The strongest correlation is expected to be 
found when Pi is constant for different pulsars. Such a corre- 
lation would be a very important observational restriction on 
pulsar emission models. 

A significant correlation between P3 and t he pulsar age has 
been reported in the past dWolszczanlll980t lAshworfhlll982l 
lRankinll9 86). As one can see in Fig.ll2lthere is no clear corre- 
lation present in our data, which is confirmed by y 2 -fitting. Th e 
figure looks qualitatively the same as Fig. 4of lRankrnl ([l986). 
although she only plots "conal P3 values". As discussed in Sect. 
14.31 it is not clear if one should make this distinction. 




t (in years) 

Fig. 12. The measured value of the vertical drift band sepa- 
ration P3 versus the pulsar age of all the pulsars showing the 
drifting phenomenon. The coherent drifters are the filled cir- 
cles, the Dif drifters (with drift feature clearly separated from 
the alias borders) are the open circles and the Dif* drifters are 
the small dots. 



There is no correlation found between Pj and the magnetic 
field s tren gth, contrary to w hat was reported by IWolszczar] 
dl980l) and lAshworfh1 (ll982). Also there seems to be no cor- 
relation between Pj, and the pulse period (consiste nt with 
Wols zczanll 19801 contrary to the tendency reported by Backer 
119731) . The evidence for a pulsar subpopulation located close to 
the P 3 = 2P Nyquist limit dWrightll20oi lRankinl H986) also 
seems to be weak. 

In a sparking gap model one would expect that the spark- 
associated plasma columns drift because of an E x B drift, 
which depends on b oth the pulse peri od and its derivative (e.g. 
iRuderman & Sutherland! fl975t Idil & Sendvkl l2000). The ab- 
sence of any correlation between P3 and a physical pulsar pa- 
rameter is difficult to explain in this model, unless many pulsars 
in our sample are aliased. If a pulsar is aliased a higher E x B 
drift can result in a lower P3 value and visa versa, making P3 
not a direct measure of the E x B drift. In the plot a distinction 
is therefore made between the coherent drifters, the diffuse Dif 
drifters and the Dif* drifters (the latter are probably more likely 
to be aliased), but this separation does not reveal a significant 
correlation. Also if Po is highly variable from pulsars to pulsar, 
any correlation with P3 is expected to be weaker. 

There are a number of pulsars found that show evidence 
for drift reversals. This kind of behavior is known for PSR 
B23 03+30 which s hows drift reversals around P3 = 2Pq 
(Red man et al.l2005l) . This is confirmed in the 2DFS of our ob- 
servation and two other other pulsars that very clearly show the 
same kind of behavior in its 2DFS are found: PSR B2148+63 
and PSR B23 10+42. Anothe r pulsar that is known to show drift 
revers als is PSR B0826-34 dGupta et alJl2004l Esamdin et alJ 
2005). This pulsar continuously changes the apparent drift di- 
rection via longitude stationary subpulse modulation. There are 
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a few more pulsars found which possibly show the same kind of 
subpulse behavior: PSRs B0037+56, B1642-03, B1944+17, 
B21 10+27, B2351+61. The evidence for drift reversals is less 
clear than for the afore mentioned P3 =s 2Pq pulsars. If pul- 
sars are proven to show drift reversals via longitude stationary 
subpulse modulation and one believes that the physical drift 
direction of the emission entities cannot change drift direction, 
it would imply that the the drifting in both drift directions is 
aliased. 

A correlation between the drift direction and the pulsar 
spin-down was reported bv lRitchings & LvneM 1975b . such that 
a high spin-down is correlated with positiv e dr ifting. The sam e 
trend was also found bv lAshworthl Jl982l) and lBackusI {l981), 
although its significance was less. In our sample there is no 
significant correlation between the drift direction and the pul- 
sar spin-down. Also there is no significant difference between 
the number of positive and negative drifters. 

The value of Pj is a much better defined parameter than 
Pi if the drift bands are non linear. Although this makes it dif- 
ficult to find any correlations with P%, the fact that the drift 
bands can be non linear is very interestin g by itself. The drift 
bands of for instance PSR B0809+74 JEdwards & Steppers! 
l2003bHProszvnski & Wolszczanlll986llWolszczan et alJll98lh 
and PS R B0320+39 JEdwards et alJ2003HEdwards & Steppers! 
2003b) show subpulse phase steps and we find that the new 
drifter PSR B2255+58 also shows a phase step. Although the 
carousel model predicts curved drift bands under certain ge- 
ometries, it cannot explain sharp discontinuities of this type. A 
way out could be that the observed discontinui ties are ca used 
by lo cal irregularities of the magnetic field ( Wols zczan et alJ 
Il98!h . 

It is argued bv lEdwardsetail (2003) that the local reduc- 
tion of the modulation index accompanied by a rapid swing 
in the modulation phase profile are the result of interference 
between two superposed drifting subpulse signals that are out 
of phase. In the non-radial pulsation s mode l subpulse phase 
steps could be explained llClemens & Rosenll2004l) . but it has 
trouble explaining how the modulation phase profile can be 
anything but piecewise linear. Visually curved drift bands are 
expected by applying subpulse amplitude windowing (pulse 
longitude dependent subpulse intensities), as pointed out by 
Iciemens & Rosenl2004 However it is observed that the phase 
profile of, for example, PSR B0818-13 makes a swing in the 
middle of the pulse profile. The phase profile is not sensitive 
for subpulse amplitude windowing and therefore the subpulse 
phase swing cannot be explained by the non-radial pulsations 
model. 

5. Summary and conclusions 

Subpulse modulation is shown to be very common for radio 
pulsars. Of the 187 analyzed pulsars 170 are shown to exhibit 
subpulse modulation. The measured upper limits on the mod- 
ulation index of many of the pulsars that do not have a mea- 
sured modulation index were high, indicating that pulsars with- 
out any subpulse modulation are probably rare. The number of 
pulsars that are known to show the drifting phenomenon is sig- 
nificantly expanded by 42. Our sample of pulsars is not biased 



on pulsar type or any particular pulsar characteristics, which al- 
lows us to do meaningful statistics on the drifting phenomenon. 

As the drifting phenomenon is thought to be exclusively a 
conal phenomenon, the modulation index of the drifters is ex- 
pected to be on average higher than the modulation index of 
the nondrifters. The absence of such a correlation (and pos- 
sible opposite trend for the coherent drifters) seems to sug- 
gest that drifting is not exclusively related to conal emission. 
Furthermore a number of pulsars classified as core single stars 
are found to show drifting, which stresses the importance to 
be unbiased on pulsar type when studying the drifting phe- 
nomenon. 

Of the 187 analyzed pulsars 68 are shown to exhibit the 
drifting phenomenon (of which 30 drift coherently or have 
drift features clearly separated from the alias borders), which 
means that at least one in three pulsars show drifting. If the 
observations that had little chance of detecting the drifting phe- 
nomenon because of an insufficient S/N are ignored, it is shown 
that at least some 55% of the pulsars drift. This implies that 
the physical conditions required for the drifting mechanism to 
work cannot be very different than the required physical con- 
ditions for the emission mechanism of radio pulsars, which is 
consistent with the absence of a strong correlation between the 
drifting phenomenon and the magnetic field strength. It could 
well be that the drifting phenomenon is an intrinsic property of 
the emission mechanism, although drifting could in some cases 
be very difficult or even impossible to detect. 

The set of modulation indices of our sample of pulsars is 
not shown to be inconsistent with four complexity parame- 
ters as derived for different emission models. Therefore none 
of the models can be ruled out based on the present observa- 
tions. Other correlations are found which should be explained 
by emission models. The population of pulsars that show the 
drifting phenomenon are on average older than the population 
of pulsars that do not show drifting and it seems that drifting is 
more coherent for older pulsars. 

Although significant correlations between P3 and the pul- 
sar age, the magnetic field strength and the pulse period have 
been reported previously, we find no such correlations in our 
enlarged sample. In a sparking gap model one would expect 
that the subpulses drift because of an E x B drift, which de- 
pends on both the pulse period and its derivative. The absence 
of a correlation between P3 and any physical pulsar parameter 
is difficult to explain in such a model, unless many pulsars in 
our sample are aliased or if P2 is highly variable from pulsar to 
pulsar. 

No significant correlation is found between the modulation 
index and the pulsar age. This seems to suggest that a high pul- 
sar age and a low modulation index are two independent factors 
for pulsars that affect the likelihood of them exhibiting coher- 
ently drifting subpulses. The evolutionary trend found seems 
to suggest that the mechanism that generates the drifting sub- 
pulses gets more and more stable as the pulsar ages, which 
could be because the magnetic axis and the rotation axis be- 
comes more aligned for older pulsars. 

The presence of subpulse phase steps results in a minimum 
in the longitude resolved modulation index. If subpulse phase 
steps are exclusively (or at least more likely) to occur in pul- 
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sars with coherently drifting subpulses, the modulation index 
of coherent drifters is expected to be on average lower. This is 
indeed the trend the we observe. It is argued bv lEdwards et alJ 
(2003) that the local reduction of the modulation index accom- 
panied by a rapid swing in the phase angle are the result of 
interference between two superposed drifting subpulse signals 
that are out of phase. It is not unlikely that interference can 
only occur if the drifting is coherent. Many pulsars are shown 
to have a modulation index which is significantly lower than 
what is expected for a purely sinusoidal subpulse signal. This 
implies the presence of a non-varying component in the pul- 
sar signal, which could be caused by superposed out of phase 
subpulse signals. 

Another possible scenario to explain the trend is that co- 
herent drifting indicates that the electrodynamical conditions in 
the sparking gap are stable, which could cause the subpulses to 
have more equal intensities. Another explanation for this trend 
would be that refraction in the magnetosphere is stronger for 
pulsars that do not show the drifting phenomenon coherently. 
In that scenario the organized drifting subpulses are refracted 
in the magnetosphere, causing the coherent drifting to be dis- 
torted. Furthermore it is expected that refraction would cause 
the subpulses to appear more disordered in the pulse window 
and that the intensity distribution of the subpulses becomes 
broadened because of lensing and possible focusing of the ra- 
dio emission. 
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Table 2. The details of all the analysed pulsars. The classification of the pulsar in the second column, where "Coh" is a coherent drifter, "Dif" 
and "Dif" are diffuse drifters with or without drift features which are clearly separated from the alias borders and "Lon" are pulsars showing 
longitude stationary subpulse modulation. The next columns are the pulse period, its dimensionless time derivative, the number of pulses in the 
observation, the signal to noise ratio, the minimum in the longitude resolved modulation index, the minimum detectable modulation index, the 
horizontal and vertical driftband separation and the figure number. 
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0.41 






IaTI 


B1811+40 




0.9311 


2.5 


10" 


15 


7200 


109 


0.56±0.1 


0.53 






ro 


J1812-2102 




1.2234 


2.4 


10" 


14 


1435 


63 




0.65 






IaTI 


B1813-17 




0.7823 


7.3 


10" 


15 


1096 


50 


0.7±0.1 


0.58 






IaTI 


B1815-14 




0.2915 


2.0 


10" 


15 


2579 


175 


0.27±0.03 


0.18 






IaTI 


B1817-13 




0.9215 


4.5 


10" 


15 


1031 


46 


0.65±0.06 


0.41 






IaTI 


B1818-04 




0.5981 


6.3 


10" 


15 


7681 


353 


0.40±0.04 


0.16 






IaTI 


B1819-22 


Dif 


1.8743 


1.4 


10" 


15 


1096 


258 


0.76±0.01 


0.10 


-9.1 H 

-1.0 


9.8+0.6 


IaTI 


B1820-11 




0.2798 


1.4 


10" 


15 


3021 


97 


0.61±0.04 


0.35 




IaTI 


B1821+05 




0.7529 


2.3 


10" 


16 


1080 


24 




0.99 






lA8l 


B1821-11 




0.4358 


3.6 


10" 


15 


1568 


59 




0.37 






lA8l 


B1821-19 




0.1893 


5.2 


10" 


15 


4105 


222 


0.632±0.004 


0.08 






lA8l 


B 1822-09 


Dif 


0.7690 


5.2 


10" 


14 


1180 


409 


0.672±0.005 


0.24 


50+!i° 


8+4 


IA.21I 


B 1822- 14 




0.2792 


2.3 


10" 


14 


3061 


79 


1.8±0.2 


0.48 




lA8l 


B1826-17 




0.3071 


5.6 


10" 


15 


2764 


361 


0.59±0.01 


0.11 






lA8l 


J1828-1101 




0.0721 


1.5 


10" 


14 


24167 


64 




0.98 






lA8l 


B 1829-08 




0.6473 


6.4 


10" 


14 


1311 


132 


0.56±0.09 


0.22 






lA8l 


B 1830-08 




0.0853 


9.2 


10" 


15 


9412 


124 


1.4+0.2 


0.22 






lA8l 


J1830-1135 


Dif 


6.2216 


4.8 


10" 


14 


1004 


140 


1.1+0.1 


0.21 


25 + _f 


2.1+0.2 


lA8l 
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Pulsar 


Class 




P 


Pulses 


S/N 


m 


^thresh 


P 2 (deg) 


Pi (Po) 


Figure 
























B1831-03 




0.6867 


4.2 


10" 




1236 


88 


0.45±0.06 


0.26 






"A9 


B1831-04 




0.2901 


7.2 


10" 




11255 


114 


0.41±0.05 


0.36 






roi 


B 1832-06 




0.3058 


4.0 


10" 




4557 


33 




0.98 








B 1834-04 




0.3542 


1.7 


10" 




4114 


73 


0.57±0.05 


0.41 








B1834-10 




0.5627 


1.2 


10" 




1286 


151 


0.32±0.01 


0.12 








J1835-1106 




0.1659 


2.1 


10" 




5125 


108 


0.89±0.04 


0.42 








B 1839+09 




0.3813 


1.1 


10" 




2056 


79 


0.48±0.06 


0.41 








B1839+56 




1.6529 


1.5 


10" 




1031 


53 


1.07±0.06 


0.44 








B 1839-04 


Coh 


1.8399 


5.1 


10" 




1025 


179 


0.49±0.04 


0.15 


-35 ** 

-15 


12.4±0.3 
























120 + ?° 


12.4±0.3 




J1839-0643 




0.4495 


3.6 


10" 




3590 


50 


0.51±0.08 


0.47 






B1841-04 


Coh 


0.9910 


3.9 


10" 




1551 


69 


0.62±0.09 


0.44 


-13 +n 

-5 


8.5±0.1 


roi 


B1841-05 




0.2557 


9.7 


10" 




3077 


115 


0.26±0.04 


0.26 








B 1842+ 14 




0.3755 


1.9 


10" 




3073 


77 


0.67±0.06 


0.46 








B 1842-04 




0.4868 


1.1 


10" 




1556 


68 


0.53±0.07 


0.45 






IaHtji 


B 1844-04 


Coh 


0.5978 


5.2 


10" 




1240 


115 


0.38±0.03 


0.26 


80 + ™ 

-45 


12±1 


lAJOl 


B 1845-01 


Dif 


0.6594 


5.3 


10" 




1031 


213 


0.28±0.01 


0.11 


90 + ™ 

-35 


14±2 




J1845-0743 




0.1047 


3.7 


10" 




8193 


115 


1.02±0.09 


0.26 






B 1846-06 


Lon 


1.4513 


4.6 


10" 




1208 


128 


1.31+0.03 


0.24 


? 


3±1 


IA.10I 


B1848+13 




0.3456 


1.5 


10" 




1544 


42 


0.55±0.09 


0.55 








B 1849+00 




2.1802 


9.7 


10" 




1030 


94 


0.25±0.02 


0.17 








Jl 850+0026 




1.0818 


3.6 


10" 




1029 


41 


0.61±0.08 


0.38 






COol 


Jl 852+0305 




1.3261 


1.0 


10" 




1323 


10 




2.21 






IA.11I 


J1852-2610 




0.3363 


8.8 


10" 




2526 


42 




0.50 






[on 


B1855+02 




0.4158 


4.0 


10" 




2051 


77 


0.7±0.1 


0.49 






IA.11I 


B1855+09 




0.0054 


1.8 


10" 


20 


158193 


175 


0.6±0.1 


0.46 






ia~tti 


B1857-26 


Dif 


0.6122 


2.0 


10" 




1339 


360 


0.47±0.02 


0.20 


803° 


7.5±0.2 


IA.21I 






















180+2? 

-45 


7.3+0.2 




B1859+03 




0.6555 


7.5 


10" 


5 


1552 


182 


0.35±0.05 


0.12 




IaTTI 


B1859+07 




0.6440 


2.3 


10" 




1030 


46 


0.47±0.08 


0.52 






IaTTI 


B 1900+01 


Dif 


0.7293 


4.0 


10" 




1171 


327 


0.56±0.05 


0.09 


30 i° 

-3 


3.4±0.7 


IaTTI 


B 1900+05 




0.7466 


1.3 


10" 




1025 


54 


0.49±0.06 


0.40 




ia~tti 


J1901-0906 


Coh 


1.7819 


1.6 


10" 




785 


147 


0.47±0.08 


0.13 


-8 

-5.6 + »i 

-0.7 


6.9±0.3 
3.1+0.1 


IA.21I 


B 1903+07 




0.6480 


4.9 


10" 




2246 


32 




0.72 






ia~tti 


B 1905+39 




1.2358 


5.4 


10" 




1032 


50 


0.78±0.07 


0.57 






IaTTI 


B 1907+00 




1.0169 


5.5 


10" 




1030 


40 


0.63±0.06 


0.44 






ia~tti 


B1907+10 




0.2836 


2.6 


10" 


5 


3023 


95 


0.51±0.07 


0.44 






lA~12l 


B1910+20 




2.2330 


1.0 


10" 




1025 


43 


0.67±0.07 


0.39 






IA.12I 


B1911+13 




0.5215 


8.0 


10" 




2055 


56 


0.32±0.04 


0.26 






lAl2l 


B1911-04 


Dif 


0.8259 


4.1 


10" 




2082 


616 


0.287±0.004 


0.08 


70 is 

-30 


15+5 


lA7l2l 


B1914+09 




0.2703 


2.5 


10" 




6145 


97 


0.59±0.08 


0.52 




IA.12I 


B1914+13 




0.2818 


3.6 


10" 


15 


4096 


185 


0.41±0.02 


0.19 






Ia7T21 


B1915+13 




0.1946 


7.2 


10" 




1541 


41 




0.46 






lA~12l 


B1916+14 




1.1810 


2.1 


10" 




1226 


23 




0.79 






IA.12I 


B1917+00 


Dif 


1.2723 


7.7 


10" 




4106 


135 


0.69±0.08 


0.27 


70 + ^ 5 

-25 


7.8±0.2 


lAl2l 


B1918+19 




0.8210 


9.0 


10" 


16 


1550 


60 


0.58±0.07 


0.48 






lAl2l 


B1919+21 


Dif 


1.3373 


1.3 


10" 




1033 


207 


0.30±0.04 


0.18 


_^ 4+0-3 
"'•^-0.2 

-11 + ! 

-1 


4.4+0.1 
4.4±0.1 


IAT221 


B1920+21 




1.0779 


8.2 


10" 




1025 


74 


0.40±0.06 


0.30 






lA~13l 


B 1924+ 16 




0.5798 


1.8 


10" 




2514 


87 


0.76±0.06 


0.44 






IA.13I 


B1929+10 


Dif 


0.2265 


1.2 


10" 




1324 


681 


0.466±0.002 


0.04 


90 +i4o 
-160 + !° 

-100 


9.8+0.8 
4.4+0.1 


IAJ31 


B1933+16 


Dif 


0.3587 


6.0 


10" 




2613 


2081 


0.240+0.001 


0.02 


300 + ™ 

-50 


2.4+0.3 


lA~13l 


B1935+25 




0.2010 




1U 


16 


7173 


66 


0.58±0.06 


0.47 




lA~13l 


B1937+21 




0.0016 


l.i 


10" 


19 


600160 


339 




0.19 






lA~13l 


B 1937-26 


Lon 


0.4029 


9.6 


10" 


16 


2129 


91 


0.93±0.08 


0.31 


? 


2.5±0.5 


IaHjI 


B 1943-29 




0.9594 


1.5 


10" 


15 


1319 


43 


0.8±0.1 


0.56 






IAJ31 


B 1944+ 17 


Dif 


0.4406 


2.4 


10" 


17 


2056 


278 


1.19±0.02 


0.17 


-11 7 +03 

11 -'-0.3 


13±5 


lAl3l 


B 1946+35 


Lon 


0.7173 


7.1 


10" 


15 


1161 


176 


0.43±0.01 


0.12 


? 


33±2 


IAT221 
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Pulsar 


Class 




P 


Pulses 


S/N 


m 


^thresh 


P 2 (deg) 


Pi (P ) 


Figure 
























B1952+29 


Dif 


0.4267 


1.7 


10" 




1536 


312 


0.46±0.06 


0.07 


-190 + !? n 

-170 


13.1+0.2 


IA.22I 






















-40 +* 

-30 


12.5+0.4 




B1953+50 


Dif 


0.5189 


1.4 


10" 




1635 


464 


0.9±0.2 


0.08 


-0.2 


25+7 


roil 


B2000+32 




0.6968 


1.1 


10" 




2060 


90 


0.8+0.1 


0.41 






ro4i 


B2000+40 


Coh 


0.9051 


1.7 


10" 




768 


109 


0.53+0.04 


0.28 


-7.3 

-0.6 


2.5+0.1 


ro2i 






















-17 + l 

-8 


2.5+0.1 




B2002+31 




2.1113 


7.5 


10" 




1115 


214 


0.54+0.02 


0.12 






[Q4l 


B2003-08 




0.5809 


4.6 


10" 




1285 


37 


0.8+0.2 


0.74 








B2011+38 


Lon 


0.2302 


8.9 


10" 




3722 


271 


1.34+0.02 


0.21 


? 


30+15 


IaT41 


B2016+28 


Dif 


0.5580 


1.5 


10" 




25899 


511 


0.61+0.01 


0.15 


' -25 


20+2 


IA.22I 






















-12+1 

-8 


4.0+0.4 




B2020+28 


Dif 


0.3434 


1.9 


10" 




3759 


1438 


0.15+0.02 


0.04 


-55 + f, 

-15 


2.5+0.2 


IX2H 






















25 + > 5 


2.3+0.1 




B2021+51 


Dif* 


0.5292 


3.1 


10" 




20326 


1180 


0.50+0.03 


0.08 


10 + ! 5 

-1 


5.1+0.3 


lAJ4l 


B2022+50 




0.3726 


2.5 


10" 




2278 


96 


0.52+0.04 


0.31 






lAJ4l 


B2043-04 


Coh 


1.5469 


1.5 


10" 




1025 


122 


0.69+0.02 


0.23 


4.5 + L 

-0.3 


2.74+0.05 


IA.14I 


B2044+15 


Dif 


1.1383 


1.8 


10" 




768 


51 


0.72+0.06 


0.40 


-7. i + ?i 

-1.4 


18+6 


lAT4l 


B2045-16 


Coh 


1.9616 


1.1 


10" 




384 


44 


0.8+0.1 


0.55 


1? +18 


3.2+0.1 


IaT41 


B2053+36 


Dif 


0.2215 


3.7 


10" 




3594 


173 


0.59+0.02 


0.22 


ioo!^° 


11+2 


IA.14I 


B2106+44 


Lon 


0.4149 


8.6 


10" 




2067 


292 


0.52+0.04 


0.19 


? 


19+7 


lA~T5l 


B21 10+27 


Dif 


1.2029 


2.6 


10" 




1032 


370 


1.06+0.02 


0.09 


140 + !* 

-15 


4.4+0.1 


Ia~T31 


B21 11+46 


Dif 


1.0147 


7.1 


10" 




3083 


258 


0.62+0.05 


0.17 


-180!* 


2.7+0.4 


IA.23I 


J2 145-0750 




0.0161 


3.0 


10" 


20 


224536 


190 


0.53+0.07 


0.16 




lAT5l 


B2148+52 




0.3322 


1.0 


10" 




2055 


89 


0.39+0.06 


0.34 






|A~T5l 


B2148+63 


Dif 


0.3801 


1.7 


10" 




2256 


143 


0.89+0.02 


0.26 


-35 + ?« 


2.4+0.3 


IaH31 


B2154+40 


Dif 


1.5253 


3.4 


10" 




1148 


475 


0.60+0.01 


0.11 




3.1+0.8 


lA"23l 


B2217+47 




0.5385 


2.8 


10" 


15 


1592 


191 


0.52+0.03 


0.18 




Ia~!31 


B2224+65 




0.6825 


9.7 


10" 




2548 


56 


0.96+0.07 


0.42 






IA.15I 


B2255+58 


Dif 


0.3682 


5.8 


10" 




2305 


506 


0.52+0.01 


0.12 




10+1 


lAT5l 


B2303+30 


Coh 


1.5759 


2.9 


10" 




1109 


96 


0.72+0.08 


0.26 


15.0+L 

-0.3 


2.1+0.1 


IaH31 


B2306+55 




0.4751 


2.0 


10" 


16 


3596 


45 


0.8+0.1 


0.70 




lA~T5l 


B23 10+42 


Coh 


0.3494 


1.1 


10" 




5006 


2073 


0.46+0.01 


0.11 


60 + ?° 

^ - 1 


2.1+0.1 


IA.23I 






















6.0 + " 

-0.3 


2.1+0.1 




B2315+21 




1.4447 


1.0 


10" 




1024 


44 


0.7+0.1 


0.58 




lAT6l 


B2319+60 


Coh 


2.2565 


7.0 


10" 




1041 


1652 


0.923+0.004 


0.03 


70 + ^ 

-10 


7.7+0.4 


lA^23l 






















12 + f 

-l 


7.7+0.1 




B2323+63 




1.4363 


Z.o 


i n~ 

1U 


15 


1036 


29 


1.6+0.1 


0.67 






lAJ6l 


B2324+60 


Dif 


0.2337 


3.5 


10" 


16 


3633 


289 


0.95+0.01 


0.16 


160!* 


17+4 


lA~T6l 


B2327-20 




1.6436 


4.6 


10" 


15 


1024 


125 


0.55+0.04 


0.15 




lAJ6l 


B2334+61 




0.4953 


1.9 


10" 


13 


3085 


19 


1.3+0.2 


1.07 






lA"T6l 


J2346-0609 


Dif 


1.1815 


1.4 


10" 


15 


770 


79 


0.77+0.04 


0.22 


-90!* 


2.3+0.1 


lA^4l 


B2351+61 


Dif 


0.9448 


1.6 


10" 


14 


7655 


661 


1.07+0.01 


0.07 


60 !f 


17+6 


lAJ6l 
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